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FOREWORD 

This  report  is  published  by  The  Aerospace  Corporation.  El  Segundo, 
California,  under  Air  Force  Contrac  t  No.  F04701  -73-C  -0074.  The  document 
was  compiled  and  edited  by  J.  A.  Pearson;  the  authors  include  R.  W.  Bruce, 

W.  A.  Feess,  E.  G.  Fotou,  J.  A.  Pearson,  and  K.  E.  Young,  Engineering 
Science  Operations,  and  V.  T.  Chiu  and  A.  B.  Prag,  Laboratory  Operations. 

This  report  documents  experimental  work  carried  out  from  June  1966 
through  June  1970.  The  present  compilation,  published  in  the  interests  of 
wider  dissemination  of  this  technical  infor mation,  was  completed  in  March 
1973  and  submitted  on  1  July  1973  to  Lt.  E.  L.  Lockwood.  SAMSO  (DYAE) 
for  review  and  approval. 

This  work  gathers  together  under  one  cover  a  complete  description 
of  the  LOGACS  orbital  accelerometer  experiment  and  the  results  obtained 
from  that  experiment.  The  editor  gratefully  acknowledges  the  assistance 
and  support  of  members  of  both  the  Upper  Atmosphere  Physic  s  Department 
of  the  Space  Physics  Laboratory,  Laboratory  Operations,  and  the  Satellite 
Navigation  Department  of  the  Electronic  s  and  Optics  Division,  Engineering 
Science  Operations,  in  the  preparation  of  this  document.  The  editor  also 
thanks  the  authors  of  the  individual  papers  included  in  these  volumes  for  per¬ 
mitting  their  work  to  be  included  and  for  their  valuable  technical  assistance. 

A  list  of  the  LOGACS  documents  is  included  in  the  bibliography  of  this  volume. 
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ABSTRACT 


The  LOGACS  (LOw-G  Ac  celerometer  Calibration  System)  experiment, 
which  contained  a  miniature  electrostatic  accelerometer  (MESA),  was  placed 
in  a  low-altitude  polar  orbit  on  22  May  1967.  The  experiment  provided 
approximately  100  hr  of  acceleration  data  from  which  the  accelerometer 
scale  factor,  accelerometer  null  bias,  and  atmospheric  drag  on  the  satellite 
were  calculated. 

Many  reports  were  written  during  1967  and  1968  on  the  LOGACS 
experiment  and  on  the  analyses  that  were  performed  on  the  data.  These 
reports  have  been  combined  into  this  document,  which  will  be  the  definitive 
report  on  the  LOGACS  experiment  by  The  Aerospace  Corporation. 

The  mission  description,  the  characteristics  of  the  accelerometer, 
the  instrumentation  used  in  the  experiment,  and  a  summary  of  the  results  of 
the  inflight  accelerometer  calibration  are  presented  in  Volume  I.  The 
accelerometer's  sensing  environment  is  discussed,  as  are  the  data  reduction 
methods  used  to  extract  the  scale  factor  and  bias  of  the  accelerometer  and 
the  methods  of  determining  the  acceleration  due  to  aerodynamic  drag. 

Volume  II  contains  data  plots  of  the  reduced  accelerometer  time 
histories  and  of  the  deduced  atmospheric  densities  over  the  flight.  Peculiar¬ 
ities  of  the  density  data  are  examined  in  some  detail  and  extensive  comparisons 
of  the  data  to  vari''1  ■;  model  atmospheres  are  made.  These  analyses  are  the 
more  interesting  because  of  the  intense  solar  and  geomagnetic  activity  during 
the  latter  part  of  the  LOGACS  flight,  which  enormously  affected  the  atmosphere. 

Models  of  atmospheric  density  have  been  developed  over  the  altitudes 
for  which  LOGACS  data  were  available.  Both  midlatitude  models  and  polar 
and  auroral  models  are  formulated  and  presented. 

An  analysis  was  also  performed  on  the  LOGACS  data  to  determine 
the  wind  magnitude  and  characteristics  normal  to  the  orbit  olane  of  the  experi¬ 
ment.  The  results  confirm  that  the  earth's  atmosphere  (up  to  the  altitude  of 
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100  nnn)  have  high-velocity  winds  that  are  present  in  high-latitude  regions 
as  a  result  of  the  great  magnetic  storm.  An  i,».;lysis  of  these  wind  data 
shows  they  can  be  described  as  convective  motion  due  to  excessive  heating 
of  the  polar  thermosphere. 
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SEC  TION  I 


LOGACS  EXPERIMENT 

A.  HISTORICAL  BACKGROUND 

It  has  long  been  recognized  that  accurate  acceleration  measurements 
aboard  satellite  vehicles  in  near-earth  orbits  would  improve  on-orbit  naviga¬ 
tion.,  would  allow  the  determination  of  a  more  accurate  atmosphere  model, 
and  would  enhance  orbit  prediction  capability.  Unfortunately,  it  has  also  been 
recognized  thav  the  accelerations  that  must  be  sensed  are  generally  several 
orders  of  magnitude  below  those  which  can  be  measured  accurately  on  earth. 
As  a  result,  even  though  an  accelerometer  might  be  capable  of  sensing  the 
extremely  low  accelerations  encountered  aboard  a  satellite  vehicle,  the 
calibration  of  the  accelerometer  prior  to  launch  could  not  be  accomplished 
except  at  levels  several  orders  of  magnitude  above  those  which  would  be 
experienced  on-orbit.  This  impass'*  led  directly  to  the  formulation  of  a  plan 
to  calibrate  an  extremely  sensitive  accelerometer  in  the  orbital  environment. 
The  result  of  this  plan  was  the  LOGACS  (LOw-G  Accelerometer  Calibration 
System)  experiment. 

Briefly,  the  plan  called  for  an  orbital  experiment  wherein  it  would 
be  demonstrated  that  the  bias  and  scale  factor  coefficients  i  f  the  accelerom¬ 
eter  could  be  determined  while  on-orbit  and,  further,  that  the  accelerometer 
would  measure  (along  the  accelerometer's  sensitive  axis)  the  nonconservative 
forces  acting  on  the  satellite  vehicle.  Measurements  would  be  made  nearly 
continuously  while  on-orbit. 

Selection  of  a  suitable  accelerometer  did  not  prove  to  be  a  dulicult 
task.  Schedule  and  limited  funding  demanded  that  the  instrument  already  be 
developed.  A  brief  survey  revealed  that  a  device  that  promised  the  required 
sensitivity  had  already  been  developed  and  produced  in  prototype  quantity  by 
Bell  Aerosystems  Company. 


Bell  Aerosystems  Company  had  initiated  in  1959,  under  Contract 
AF  33(616) -6637,  a  program  to  develop  a  single  axis,  force  rebalance, 
digital,  electrostatic  accelerometer  (ESA)  for  extremely  low  acceleration 
measurements.  The  force  rebalance  was  to  be  accomplished  by  sensing  the 
displacement  of  a  proof  mass  and  using  that  displacement  signal  to  close  the 
loop  to  an  electrostatic  force  generator  to  null  the  displacement  signal.  The 
feasibility  of  the  design  was  demonstrated  after  four  years  of  development 
effort.  Under  sponsorship  of  the  NASA  Lewis  Research  Center,  a  contract 
(NAS-3 -4102)  was  initiated  in  196  3  to  develop  two  miniature  versions  of  the 
electrostatic  accelerometer  (MESA  IA  and  MESA  IB).  The  MESA  LA  was 
successfully  qualified  to  the  Thor -Delta  launch  environment.  The  MESA  IB 
was  a  low  power  version  which  was  successfully  qualified  to  the  Atlas-Agena 
launch  environment. 

Since  it  was  recognized  that  a  ride  as  a  noninterfering  passenger 
aboard  an  orbiting  Agena  vehicle  could  be  arranged,  it  was  decided  to  pro¬ 
ceed  with  the  plan.  The  preliminary  design  of  the  experiment  was  performed 
by  Ihe  Aerospace  Corporation.  Discussions  were  held  with  Lockheed 
Missiles  and  Space  Company  (LMSC)  and  Bell  Aerosystems  Company,  and 
subsequently,  LMSC  was  authorized  to  proceed  with  the  acquisition  of  the 
necessary  equipment  to  fly  two  LOGACS  experiments  aboard  separate  Agena 
vehicles  to  be  designated  by  the  Space  Systems  Division  of  the  USAF. 

B.  EXPERIMENT  DESIGN 

The  fundamental  experiment  design  problem  was  the  provision  of  a 
means  to  calibrate  the  accelerometer  bias  and  scale  factor  while  on-orbit. 

Let  us  define  bias  and  scale  factor  as  follows 

A  =  K(B  +  A  ) 
m  s 

See  Section  C  for  a  detailed  description  of  the  instrument  and  its  operation. 
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whe  re 


=  accelerometer  output  (pulses/' sec) 

K  =  scale  factor  (pulses/sec-g) 

B  =  bias  (g) 

=  sensed  acceleration  (g) 

The  calibration  of  accelerometer  bias  can  best  be  illustrated  by  a 
simple  example.  Let  us  fi:’  the  accelerometer  such  that  the  sensitive  axis  is 
aligned  with  the  drag  vector  (call  this  Mode  3).  Next,  reorient  the  sensitive 
axis  by  turning  it  end  over  end  so  that  it  is  again  aligned  with  the  drag  vector 
but  with  the  opposite  sense  (call  this  Mode  4).  Let  us  further  assume  that 
th'1  drag  vector  is  constant  during  the  measurement  interval;  then  the  equations 
for  the  accelerometer  output  for  Mode  3  and  Mode  4  are 


Mode  3:  -  K(B  t  Arf  ) 


(1) 


Mode  4: 


rr>4 


K(B  -  A 


drag* 


(2) 


Adding  Eqs.  (1)  and  (2),  we  obtain 


13  = 


A  ,  +  A  . 
m3  m4 

2K 


(3) 


Thus,  if  we  know  K  (scale  factor),  we  can  solve  for  bias. 

To  solve  for  K,  W(  must  create  a  standard  linear  acceleration.  A 
centrifuge  is  a  common  solution  to  the  problem  of  creating  a  linear  accelera¬ 
tion  standard.  While  a  centrifuge  is  normally  employed  to  create  a  large 
acceleration,  there  is  no  conceptual  reason  for  noc  employing  it  to  create 
small  accelerations. 

It  was  decided  to  employ  a  small  centrifuge  which  could  be  pro¬ 
grammed  to  run  at  a  fixed  speed  (and  hence  generate  a  fixed  centripetal 
acceleration  of  magnitude  Rw^,  where  R  is  the  distance  from  the  center  of 
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mass  Of  the  sensitive  element  to  the  center  of  rotation  and  w  is  the  fixed 
angular  velocity  of  the  centrifuge).  The  centrifuge  also  could  be  pro¬ 
grammed  to  stop  at  specified  locations.  Thus  a  single  device  could  be 
employed  to  calibrate  for  bias  and  scale  factor.  The  bias  and  scale  factor 
calibration  concept  is  summarized  in  Fig.  1. 

The  manner  in  which  the  scale  factor  can  be  extracted  from  the 
accelerometer  output  data  can  be  illustrated  by  the  following.  Define  Mode  1 
as  the  condition  wherein  the  centrifuge  rotates  at  a  fixed  speed  wand  Mode  2 
as  the  condition  wherein  the  centrifuge  rotates  at  a  fixed  speed  2w.  Then 


BIAS 


ACCELEROMETER  OUTPUT  (1)  =  SCALE  FACTOR 


BIAS  +  DRAG 


BY  REVERSING  THE  SENSITIVE  AXIS 


ACCELEROMETER  OUTPUT  (2)  *  SCALE  FACTOR 

(bias  -  drag! 


SCALE  FACTOR 

MESA  MOUNTED  A  KNOWN  DISTANCE  R  FROM  TABLE 
CENTER 

TABLE  ROTATION  RATE  (w)  CONTROLLED 
Rw2  BECOMES  ACCELERATION  STANDARD 


Fig.  1.  Bias  and  Scale  Factor  Calibration  Concept 
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the  equations  which  describe  the  accelerometer  output  for  Mode  1  and 
Mode  2  are 


Mode  1  s  A 


,i  * K  [ 


B  +  A,  cos  ut  +  Rw 
drag 


Mode  2:  A 


“E 


B  +  A  ,  cos  2wt  +  4Ro 
drag 


For  the  simple  case  where  drag  is  essentially  constant,  as  at  apogee, 
we  can  imagine  integrating  the  accelerometer  output  for  one  rev.  of  the 
cer.trifuge  and  then  dividing  the  result  by  the  time  over  which  the  rev. 
occurred.  Since  A^  would  average  to  zero  over  one  rev. ,  the  average 
accelerometer  output  over  one  rev.  for  Mode  1  and  Mode  2  can  be 
expressed  as 


An,>  ■  h  C  Aml  d0  ■  K  [b  *  R-2] 

Am2  =  £  t  AmZ  d6  -  K  tB  *  4R“*2] 


If  we  subtract  Eq.  (7)  from  Eq.  (6),  we  find 


A  _  -  A  , 
m2  ml 


We  also  find,  if  we  multiply  Eq,  (6)  by  four  and  then  subtract  the  Eq.  (7) 
from  it,  that 


[4A  ,  -  A  -  "1 

— T-£nl  ™ 

^m2  ’  ^ml  J 


1-5 


-y -a-  ’iSSSSiSSXlSSf'.  WWWV"*T*'  'v--'5JWV 


It  should  be  noted  that  Eq.  (9)  does  not  provide  the  same  accuracy 
in  determination  of  instrument  bias  as  would  Eq.  (3).  This  is  because  Eq.  (9) 
is  dependent  on  a  small  difference  of  two  large  numbers  for  its  determination 
of  bias.  In  the  actual  determination  of  the  instrument  parameters,  the  drag 
is  not  constant  and  the  simplified  expressions  shown  above  do  not  apply.  The 
actual  case  is  covered  by  representing  the  drag  by  a  second-order  polynomial 
with  unknown  coefficients.  Such  a  polynomial  representation  permits  the 
development  of  equations  which  are  equivalent  to  but  more  complex  than 
Eqs.  (1)  through  (9). 

When  the  initial  calibration  concept  outlined  above  had  been  estab¬ 
lished,  it  was  necessary  to  generate  a  preliminary  design  to  the  level  of  a 
block  diagram  and  signal  flow  diagram  as  well  as  a  preliminary  estimate  of 
the  experiment  parameters.  A  schedule  and  cost  estimate  were  also  made. 

It  was  decided  to  proceed  with  the  design  and  fabrication  of  three 
LOGACS  assemblies.  One  was  to  be  used  for  qualification  testing;  the  other 
two  were  to  be  scheduled  for  flight.  The  Agena  vehicle  on  which  the  LOGACS 
would  be  mounted  would  be  stabilized  with  respect  to  the  local  vertical  and 
the  orbit  plane.  The  location  of  the  equipment  would  be  the  aft  equipment 
racks.  The  vehicle  would  fly  a  near  earth  orbit  with  a  small  eccentricity 
(i  =  0.  014). 

At  this  point,  a  contract  was  let  to  LMSC  to  acquire,  test,  integrate 
with  the  Agena,  and  support  two  LOGACS  flights.  LMSC  was  directed  to  sub  - 
contract  to  Bell  Aerosystems  Company  the  key  subsystem  of  the  LOGACS, 
which  would  consist  of  the  centrifuge  on  which  would  be  mounted  a  single 
miniature  electrostatic  accelerometer  (MESA)  and  its  electronics. 

C.  DESCRIPTION  OF  MESA 

The  basic  MESA  instrument  consists  of  an  electrostatically  supported 
proof  mass  that  is  electrostatically  pulse  rebalanced  along  a  preferred  or 
sensitive  axis.  Figure  2  show's  the  complete  MESA  instrument  while  Fig.  3 
shows  the  instrument  without  the  insulating  cover.  The  exploded  view  of  the 
MESA  shown  in  Fig.  4  and  the  cross  section  shown  in  Fig.  5  will  be  useful 
in  the  description  of  the  MESA  which  follows. 
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MESA  Without  Insulating 


The  proof  mass  consists  of  a  hollow  cylinder  with  a  centrally  located 
flange.  The  proof  mass  is  supported  by  electric  fields  generated  by  the  eight 
support  electrodes,  each  of  which  is  in  series  with  a  tuning  inductor.  The 
capacitance  in  each  series  tuned  circuit  changes  as  a  function  of  the  relative 
displacement  of  the  proof  mass  from  Lhe  support  electrode.  Since  each 
resonant  circuit  is  tuned  to  operate  on  the  high-frequency  side  of  the  reso¬ 
nance  curve,  each  electrode -to -proof-mas s  voltage  varies  directly  as  the 
electrode-to-proof-mass  displacement.  As  a  result,  the  attractive  force 
between  each  electrode  and  the  proof  mass  increases  as  the  electrode -to- 
proof-mass  distance  increases,  thus  producing  a  stable  support.  This  cross¬ 
axis  support  system  for  the  proof  mass  constrains  the  proof  mass  trans- 
lationally  along  the  two  axes  perpendicular  to  the  sensitive  axis  and 
rotationally  about  the  same  two  axes.  The  proof  mass  remains  rotationally 
unconstrained  about  the  sensitive  axis.  The  remaining  degree  of  freedom  of 
the  proof  mass  is  translation  along  the  sensitive  axis,  which  is  constrained 
by  the  pulse  rebalance  system. 

The  proof  mass  support  capability  provided  by  the  cross -axis  sup¬ 
port  system  can  be  changed  by  simply  changing  the  voltage  applied  to  the 
electrodes.  Since  the  support  capability  varies  as  the  square  of  the  electrode - 
to-pioof-mass  voltage,  it  is  apparent  that  orders  of  magnitude  change  in  the 
support  capability  can  easily  be  achieved. 

The  particular  support  capabilities  that  were  chosen  for  the  MESA 

employed  in  LOGAC5  were  (a)  1.  5  g  for  use  during  ground  test  and  on-orbit 

.  3 

by  command  and  (b)  7.5  X  10  g  for  on-orbit  use. 

The  displacement  of  the  proof  mass  along  the  sensitive  axis  is 
detected  by  an  extremely  sensitive  capacitance  bridge  which,  in  simple  form, 
is  illustrated  in  the  upper  left  corner  of  Fig.  5.  The  flange  of  the  proof  mass 
serves  as  the  moving  plate  of  the  capacitance  bridge  while  the  pickoff  ring  in 
each  housing  is  the  fixed  plate.  Each  pickoff -ring -to -flange  capacitance 
forms  a  leg  of  the  bridge  while  the  other  two  'egs  can  be  assumed  fixed  and 
equal. 
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Displacement  of  the  proof  mass  causes  unbalance  in  the  bridge,  and 
the  bridge  output,  when  amplified  and  phase  -  sensitive  demodulated,  provides 
a  signal  whose  amplitude  is  proportional  to  displacement  from  center  and 
whose  polarity  indicates  sense. 

When  the  demodulated  output  of  the  capacitance  bridge  reaches  a 
fixed  trigger  level  (E,j,  or  -E.p),  the  logic  within  the  instrument  gates  a  pre¬ 
cisely  controlled  pulse  from  the  pulse  generator  to  the  forcer  rings.  A  pulse 
always  goes  to  a  pair  (inner  and  outer)  of  forcer  rings  in  one  housing  or  the 
other. 

When  a  pulse  is  directed  to  a  forcer  ring  pair,  the  pulse  polarity  is 
such  that  a  positive  pulse  is  sent  to  one  ring  and  a  negative  (i.  e. ,  inverted) 
pulse  is  sent  to  the  other  ring.  Since  the  force  on  the  proof  mass  is  propor¬ 
tional  to  the  square  of  the  pulse  voltage,  the  force  resulting  from  the  positive 
and  negative  pulses  will  be  the  same  foi  equal  area  rings.  More  importantly, 
the  technique  of  using  opposite  polarities  for  the  pulse  on  each  forcer  ring 
prevents  the  buildup  of  charge  on  the  proof  mass. 

As  described  above,  the  simultaneous  application  of  positive  and 
negative  (i.  e.,  inverted)  pulses  to  the  inner  and  outer  forcer  rings  on  the 
same  side  of  the  proof  mass  provides  a  restoring  force  that  attracts  the 
proof  mass  to  the  forcer  rings  to  which  the  pulses  were  applied.  The  appli¬ 
cation  of  the  same  pulses  to  the  forcer  ring  pair  on  the  opposite  side  of  the 
instrument  provides  a  restoring  force  in  the  opposite  direction.  Thus, 
knowledge  of  which  pair  of  forcer  rings  is  being  pulsed  is  the  only  informa¬ 
tion  required  to  establish  direction  of  the  sensed  acceleration.  The  logic  that 
directs  the  pulses  to  one  forcer  ring  pair  or  the  other  is  set  by  the  polarity 
of  the  demodulated  capacitance  bridge  output. 

The  rate  at  which  pulses  are  demanded  is  directly  proportional  to 
the  net  force  on  the  proof  mass  and  hence  to  acceleration  applied  along  the 
sensitive  axis. 

A  more  detailed  description  and  analysis  of  the  instrument  operation 
is  given  in  Ref.  1. 
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D. 


EQUIPMENT  DESIGN 


The  longest  lead  time  items  in  the  LOGACS  equipment  were  the 

centrifuge  and  the  MESA.  Thus,  the  design  constraints  for  these  items  had 

to  be  established  first.  The  Aer*.  space  Corporation,  in  conjunction  with  the 

two  contractors,  determined  the  design  constraints  for  the  centrifuge  and  the 

MESA.  It  was  decided  that  the  MESA  would  have  a  full  scale  range  of 
-  4 

5  X  10  g,  which  would  correspond  to  an  output  of  5000  pulses/sec.  It  was 
also  decided  that  the  centrifuge  should  be  programmed  to  have  four  modes, 
each  of  512.  sec  duration.  The  four  modes  were  to  be  as  defined  in  Table  1. 
The  accelerometer  output  would  be  as  shown  in  Fig.  6. 


Table  1.  LOGACS  Centrifuge  Modes 


Mode 

Do  scription 

Accelerometer  Output 

1 

Rotation  at  fixed  angular 
speed  u- 

2 

Aj  =  K(B  +  -A-drag  cos  wt  +  ) 

2 

Rotation  at  fixed  angular 
speed  2u> 

A.,  =  K(B  +  A,  cos  2wt  +  4Rw^) 

2  drag  ' 

3 

No  rotation  (w  =  0),  sensitive 
axis  forward 

A3  =  K(B  +  Adrag) 

4 

No  rotation  (w  =  0),  sensitive 
axis  aft 

A4  =  K(B  *  Adrag» 

The  centrifuge  itself  is  shown  in  Fig.  7.  The  rotating  member 
housed  the  MESA  and  the  MESA  electronics.  The  distance  from  the  center  of 
mass  of  the  MESA  proof  mass  to  the  center  of  rotation  of  the  centrifuge  was 
3.9701  ±  0.  0024  in.  The  large  bull  gear  shown  in  Fig,  7  was  part  of  the 
rotating  member.  The  power  to  drive  the  rotating  member  was  provided  l»y 
a  hysteresis  synchronous  motor.  The  motor  was  basically  a  400-Hz  machine, 
but  it  was  driven  at  300  Hz  and  600  Hz  with  appropriate  voltages  to  provide 
the  two  speeds  for  Mode  1  and  Mode  2,  respectively.  The  power  and  signal 
leads  from  the  frame  of  the  centrifuge  assembly  to  the  rotating  member 
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MODE  3  -  4 
TRANSITION 


Characteristic  Accelerometer  Output  (constant  drag  easel 


1 XX i AC  S  Ci-ntril 


were  conveyed  through  a  slip  ring  assembly  mounted  adjacent  to  the  bull 
gear.  The  brush  block  of  the  3lip  ring  assembly  can  be  seen  in  Fig.  7.  The 
position  pickoffs  indicated  when  the  rotating  member  was  positioned  such  that 
the  MESA  sensitive  axis  was  fore  or  aft.  These  were  simply  reed  relays 
actuated  by  a  small  magnet  inset  into  the  periphery  of  the  rotating  member. 
The  reed  relays  can  be  seen  fastened  to  the  framework  of  the  centrifuge  in 
Fig.  7.  The  problem  of  lubrication  for  the  slip  rings  while  in  the  orbital 
environment  was  avoided  by  maintaining  a  partial  atmosphere  of  dry  air  in 
the  centrifuge  assembly.  The  cover  that  enclosed  the  basic  centrifuge  shown 
in  Fig.  7  contained  a  relief  valve  that  maintained  4  psi  gauge  pressure  within 
the  centrifuge  assembly.  The  seal  bet  veen  the  cover  and  the  baseplate  of  the 
centrifuge  was  effected  with  the  gasket  shown  in  Fig.  7. 

One  can  calculate  the  two  rotation  rates  (w)  and  (2w)  of  the  r  otating 
member  using  the  known  gear  ratio  of  the  gear  train  and  the  motor  excitation 
frequencies  of  300  Hz  and  600  Hz.  The  values  are 

oj  =  0.4449445  rpm 


and 


2u;  =  0.  8898890  rpm 

Since  the  experiment  clock  was  the  reference  for  the  motor  supply  fre¬ 
quencies,  and  it  is  known  that  the  experiment  clock  ran  fast  by  one  part  in 
50,  000  (see  Appendix  A),  the  above  values  should  be  multiplied  by  1.00002  to 
obtain  the  actual  values.  Performing  the  multiplication  and  rounding  off  to 
five  significant  figures  y'  Ids  the  following  values  for  (w)  and  (2w): 

u>  =  0.44495 


and 


2w  =  0.  88990 
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Use  of  the  corrected  values  rather  than  the  ideal  values  results  in  no 
iiscernible  difference  in  the  parameters  reduced  by  the  data  analysis 
p  rogram. 

As  time  progressed,  the  remainder  of  the  LOGACS  subsystems 
were  defined  and  design  or  acquisition  of  these  elements  began  to  take  place. 
The  final  block  and  signal  flow  diagram  which  evolved  is  shown  in  Fig.  B. 

It  will  be  helpful  if  a  few  explanatory  comments  are  made  about 
Fig.  8.  Note  that  the  centrifuge  (which  contains  the  MESA  and  its  electronics) 
and  the  LOGACS  I^U  (interface  unit)  have  a  clean  interface  with  the  remainder 
of  the  elements.  fhese  two  blocks  comprised  the  portion  furnished  to  LMSC 
by  Bell  Aerosystens  Company.  The  remainder  of  the  LOGACS  equipment 
was  acquired  and  integrated  with  the  Agena  vehicle  by  LMSC. 

The  primary  channel  of  information  which  was  to  be  recovered  from 
LOGACS  was  the  accelerometer  output.  As  indicated  in  Fig.  8,  the  accel¬ 
erometer  output  was  formatted  on  two  lines.  One  line  contained  the  output 
pulses  while  the  other  line  carried  a  binary  signal  that  provided  the  sense  or 
polarity  to  he  associated  with  each  output  pulse.  Thus  the  polarity  of  the 
input  acceleration  could  be  recovered  for  each  output  pulse. 

The  "pulse  line"  and  the  "polarity  line"  were  fed  to  logic  that  was 
composed  of  two  up-down  counters  and  a  register.  The  two  counters  were 
interrogated  on  alternate  1-sec  intervals.  Thus,  the  first  counter  would 
count  for  a  1-sec  interval;  then,  while  the  second  counter  was  counting  for 
the  next  succeeding  second,  the  count  magnitude  and  polarity  of  the  first 
counter  would  be  transferred  to  the  register  and  subsequently  to  the  com¬ 
mutator  and  tape  recorder.  The  1-sec  counts  were  transmitted  as  an  octal 
number.  Each  integer  of  the  octal  number  was  represented  as  a  voltage  on 
an  octal  ladder.  The  least  significant  integer  in  the  octal  number  represented 
one  accelerometer  output  pulse. 

Since  the  tape  recorder  had  two  channels  and  only  one  channel  was 
required  for  the  commutated  information,  the  second  channel  was  devottd  to 
recording  accelerometer  output  but  in  a  different  format  than  was  employed 
with  the  commutator.  The  output  pulse  line  and  the  polarity  line  were  made 
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VEHtQJE  interface 


Fig.  8.  Block  and  Signal  Flow  Diagram  for  LOGACS  Experiment 


available  to  logic  consisting  of  an  op-down  counter  that  overflowed  at  a  count 
of  eight  and  set  a  flip-flop  to  trigger  either  a  plus  pulse  or  a  negative  pulse 
depending  on  the  direction  of  the  overflow  of  the  counter.  The  same  flip-flop 
also  reset  the  counter.  This  reformatted  pulse  train  was  directly  recorded 
on  the  second  channel  of  the  tape  recorder.  The  utility  of  this  second  channel 
was  twofold:  (a)  it  offered  a  degree  of  redundancy  and  (b)  it  offered  data  with 
finer  time  resolution  than  the  1-sec  mil se  counts  that  were  placed  on  the 
commutated  channel. 

The  entire  LOGACS  eouipment  was  synchronized  from  the  airborne 
clock.  This  clock  was  a  high  accuracy  8-day  clock  with  an  octal  ladder  output 
format.  The  least  significant  octal  count  corresponded  to  4  sec.  The  incre¬ 
ment  on  the  second  octal  ladder  was  32  sec;  on  the  third,  256  sec;  on  the 
fourth,  2048  sec;  and  on  the  fifth,  16,384  sec.  A  sixth  octal  ladder  was 
reserved  for  count  days  (i.  e. ,  increments  of  86,400  sec).  Each  24  hr  after 
turn-on  of  the  clock,  each  of  the  first  five  octal  ladders  was  reset  to  zero. 

The  several  additional  functions  provided  by  the  clock  were: 

a.  Providing  a  mode  change  pulse  every  512  sec.  The  logic  for 
the  mode  change  was  related  to  even  steps  in  the  ladder  which 
had  an  increment  size  of  256  sec.  Since  each  of  the  first  five 
ladders  was  reset  to  zero  at  the  end  of  each  24-hr  period,  the 
logic  would  issue  a  mode  change  pulse  at  86,400  sec.  Since 
86,400  lacks  128  of  being  an  integral  multiple  of  512,  the  last 
mode  of  each  24-hr  period  was  not  of  512  sec  duration  but  was 
512  minus  128,  or  384  sec  long 

b.  Providing  a  2400-Hz  square  wave  as  a  reference  frequency 

c.  Providing  a  commutator  synchronizing  pulse  each  2  sec  (the 
commutator  was  a  0.  5  rps  X  60  point  electronic  commutator) 

The  clock  package  also  contained  the  l -sec  counters  and  their 
accompanying  register  as  well  as  the  up-down  eight  counter  plus  circuits  to 
format  several  miscellaneous  signals. 

The  measured  drift  of  the  airborne  clock  was  established  from  post¬ 
flight  data  analysis.  The  drift  was  found  to  be  approximately  1  sec  every 
50,  000  sec.  The  details  of  this  determination  can  be  found  in  Appendix  A. 
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The  two-channel  tape  recorder  had  approximately  163  min  of 
storage.  The  recorder  had  a  26/1  reproduce/record  speed  ratio  so  that  the 
information  acquired  during  each  orbital  rev.  could  be  transmitted  to 
the  ground  station  during  a  single  station  acquisition. 

The  LOGACS  experiment  included  its  own  telemetry  transmitters 
and  attendant  antennas  as  well  as  an  S-band  transponder  for  accurate  tracking 
of  the  vehicle.  The  LOGACS  experiment  also  included  a  complete  command 
receiving  system  so  that  commands  could  be  issued  by  the  ground  stations 
for  purposes  of  controlling  the  LOGACS  equipment  while  on-orbit. 

The  information  contained  on  the  commutator  consisted  of: 

a.  Accelerometer  pulse  counts  per  second 

b.  Experiment  clock  time  in  seconds  since  turn-on 

c.  Mode  status 

d.  Two  centrifuge  table  position  pickoff  signals 

e.  Monitor  of  Agena  attitude  control  jet  firings 

f.  Agena  gyro  output  signals 

g.  Agena  horizon  scanner  output  signals 

h.  MESA  instrument  temperature 

i.  MESA  electronics  temperature 

j.  28-V  unregulated  voltage 

k.  MESA  proof  mass  position  signal 

l.  MESA  proof  mass  suspension  signal 

m.  Necessary  synch  and  calibration  values 

The  manner  in  which  the  Agena  attitude  control  jet  firings  [item  ie) 
above]  were  recovered  is  not  obvious  when  it  is  recognized  that  the  typical 
duration  of  each  firing  is  only  0.  020  sec.  Since  the  commutator  scanned 
60  points  each  second,  some  formating  of  the  relay  closure  that  activated 
each  jet  firing  was  needed.  Each  relay  closure  was  routed  to  logic  that 
generated  a  pulse  as  shown  in  Fig.  9. 


1-20 


Fig.  9.  Relay  Closure  Logic  Pulse 


Each  of  the  sbc  jets  was  scanned  twice  on  the  commutator,  and  the 
spacing  between  the  first  and  second  placement  of  each  jet  on  the  commutator 
was  approximately  30  spaces.  Since  the  commutator  was  synchronized  with 
the  experiment  clock,  the  time  at  which  the  leading  edge  of  the  commutated 
segment  occurred  would  be  known.  Knowledge  of  the  voltage  that  existed 
during  the  scan  of  the  commutator  segment  concerned  with  a  particular  jet, 
combined  with  the  known  starting  time  of  that  particular  segment,  makes  it 
possible  to  determine  the  time  for  the  beginning  of  each  jet  firing. 

The  attitude  control  monitor  described  above  would  not  be  adequate 
if  the  attitude  control  jet  firing  rate  exceeded  1  firing  per  second  for  a  partic¬ 
ular  jet.  Postflight  analysis  has  shown  that  firing  rates  did  not  exceed  1 
firing  per  second  for  a  particular  jet  except  for  an  occasional  occurrence 
during  extremely  high  drag  regions  late  in  the  LOGACS  flight.  Further  dis¬ 
cussion  of  the  attitude  control  jet  firing  rate  is  contained  in  Section  III. 

A  list  of  the  drawings  and  documents  that  describe  the  LOGACS 
equipment  is  presented  in  Appendix  B. 

E.  PREFLIGHT  CHECKOUT 

The  entire  LOGACS  equipment  array  was  accommodated  on  the 
bottom  panel  and  the  right  side  panel  of  the  Agena  aft  equipment  rack  (see 
Fig.  10).  The  side  panel  is  shown  in  Fig.  11,  and  the  bottom  panel  is  shown 
in  '•'ig.  12.  As  can  be  seen  in  these  photographs,  reflective  tape  was  applied 
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TAPE  RECORDER 


Fig.  10.  I  ccation  of  LOGACS  Equipment  on  Agena 


Reproduced  from 
best  available  copy< 


L.OGACS  Bottom  Panel 


to  the  surfaces  of  much  of  the  equipment  to  provide  the  thermal  conditioning 
of  the  experiment  for  orbital  use.  The  antennas  and  the  antenna  ground  plane 
have  not  yet  been  installed  on  the  bottom  panel  shown  in  Fig.  12. 

After  the  panels  were  checked  out  as  an  integrated  system  at  the 
LMSC  facility  at  Sunnyvale,  California,  they  were  shipped  to  the  LMSC 
facility  at  Vandenberg  AFB.  At  Vandenberg  the  panels  were  subjected  to  a 
complete  sequence  of  tests  prior  to  integration  with  the  Agena.  The  test 
sequence  included  testing  of  the  command  equipment  by  sending  actual  com¬ 
mands  from  the  Vandenberg  launch  complex.  After  satisfactory  completion 
of  the  testing,  the  panels  were  mated  with  the  Agena  and  the  testing  sequence 
was  essentially  repeated.  At  the  conclusion  of  this  last  set  of  tests,  the 
LOGACS  equipment  was  to  remain  dormant  until  the  power  turn-on  command 
was  sent  after  injection  of  the  Agena  into  orbit.  Subsequent  to  the  final 
LOGACS  tests,  the  Agena  was  mated  to  the  Atlas  launch  vehicle  and  the 
normal  preparations  for  launch  proceeded. 

F.  FLIGHT  OPERATIONS 

After  the  Agena  was  successfully  injected  into  orbit,  power  was 
applied  to  the  LOGACS  equipment  by  means  of  a  switch  actuated  by  the  Agena 
'D"  timer. 

The  USAF  Satellite  Control  Facility  provided  the  means  whereby  the 
LOGACS  equipment  could  be  controlled  via  real-time  commands  issued  from 
the  grovind  stations  in  the  satellite  control  network.  These  same  stations  pro¬ 
vided  radar  tracking  of  the  Ager.a  on  which  the  LOGACS  equipment  resided. 

Continuous  monitoring  of  the  LOGACS  experiment  was  required, 
since  at  each  ground  station  acquisition  it  was  necessary  to  command  the 
telemetry  on  and  the  tape  recorder  to  switch  from  record  to  reproduce.  The 
telemetry  data  transmitted  from  the  LOGACS  experiment  were  recorded  on 
magnetic  tape  at  the  ground  stations.  In  addition,  the  real  -time  LOGACS  data 
transmitted  at  each  station  acquisition  were  decommutated  and  printed  out  at 
the  USAF  Satellite  Test  Center  in  Sunnyvale,  California.  Thus,  a  real-time 
monitor  of  the  status  of  the  LOGACS  experiment  was  maintained  and  some 
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preliminary  calculations  of  accelerometer  instrument  bias  and  scale  factor 
as  well  as  atmospheric  drag  were  made  during  the  100-hr  lifetime  of  the 
LOGACS  experiment. 

In  addition  to  the  normal  commands  issued  at  each  station  acquisi  ¬ 
tion,  several  special  command  situations  were  planned  during  the  LOGACS 
flight.  An  orbit  adjust  rocket  firing  was  commanded  at  the  Rev.  18.4 
station  acquisition  so  that  the  orbital  lifetime  of  the  Agena  vehicle  could  be 
extended.  At  the  Rev.  48.4  station  acquisition,  a  command  was  issued  to 
switch  to  the  higher  cross  -axis  support  voltages  for  the  MESA.  The  higher 
support  voltages  were  left  on  until  the  Rev.  56.  1  station  acquisition,  at 
which  time  a  command  was  issued  to  return  the  cross-axis  support  voltages 
to  the  normal  level  used  on-orbit.  The  purpose  of  leaving  the  higher  sup¬ 
port  voltages  on  for  several  revs,  was  to  calibrate  for  instrument  bias  and 
scale  factor  with  the  higher  support  voltages  as  well  as  the  same  quantities 
when  the  normal  support  voltages  used  on-orbit  were  employed.  The  values 
so  obtained  were  intended  to  be  used  for  comparison  with  the  ground  mea¬ 
surements  obtained  pnor  to  flight. 

G.  EQUIPMENT  PROBLEMS  DURING  THE  FLIGHT 

At  the  Rev.  1.3  station  acquisition,  the  real-time  data  were 
decommutated,  and  a  printout  was  provided  subsequent  to  the  acquisition.  It 
was  found  that  all  the  LOGACS  equipment  appeared  tc  be  operating  properly 
but  that  the  MESA  output  was  saturated;  i.  e.  ,  a  full  scale  output  was  present 
constantly.  It  had  not  been  anticipated  that  the  instrument  w^ould  remain 
saturated  for  so  long  after  launch,  but  it  was  decided  to  simply  wait  for  the 
next  acquisition  and  see  what  existed  at  that  time.  At  the  Rev.  2.4 
station  acquisition,  the  LOGACS  data  printout  was  obtained,  and  the  MESA 
output  was  again  found  to  be  at  full  scale.  Since  it  had  been  anticipated  that 
the  MESA  output  would  remain  saturated  for  just  a  few  hundred  seconds  after 
injection  into  orbit,  there  was  a  good  deal  of  anx'ety  among  the  experimenters 
at  this  time.  Though  not  planned  for  this  purpose,  the  1.5  g  support  voltage 
was  commanded  on  for  a  brief  period  at  the  next  station  acquisition  (Rev.  5.2) 
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on  the  presumption  that  if  the  test  mass  were  somehow  "stuck,"  this  would 
serve  to  free  it.  The  data  available  following  the  acquisition  showed  that 
this  action  was  unnecessary  because  the  instrument  was  operating  properly 
at  the  beginning  of  the  Rev.  5.2  acquisition.  Further  discussion  of  the  long 
saturation  period  appears  in  Section  J. 

Another  instance  of  saturation  for  a  long  period  of  time  occurred 
immediately  after  the  orbit  adjust  rockets  were  fired  at  the  Rev.  18.4 
station  acquisition.  The  orbit  adjust  rockets  consisted  of  two  small  solid 
propellant  charges  whose  net  thrust  was  directed  along  the  longitudinal  axis 
of  the  Agena.  The  result  of  firing  the  rockets  was  to  increase  the  height  of 
apogee  by  a  few  miles  and  thereby  extend  the  orbital  lifetime  of  the  Agena 
vehicle.  The  sensitive  axis  of  the  MESA  was  fixed  pointing  aft  with  respect 
to  the  Agena  during  the  orbit  adjust  rocket  firing,  and  hence  the  output  of  the 
MESA  saturated  at  this  time.  Again  the  MESA  output  remained  saturated  for 
an  extended  period  of  time.  At  the  Rev.  2C.  5  station  acquisition,  the  MESA 
output  was  found  to  be  on  scale.  No  further  instances  of  a  saturated  output 
occurred  during  the  LOGACS  flight. 

It  was  anticipated  that  the  MESA  instrument  bias  would  change  in 
accordance  with  the  ratio  of  the  support  volteges.  But  when  the  MESA  cross¬ 
axis  support  voltages  were  changed  to  the  higher  level  at  the  Rev.  48.4 
station  acquisition,  the  instrument  bias  only  changed  in  magnitude  by  about  a 
factor  of  10  compared  to  the  change  in  the  ratio  of  the  support  voltages  of  200. 
Thus,  the  bias  did  not  change  as  expected.  Further  discussion  of  this  prob¬ 
lem  will  be  left  to  Section  J. 

H.  DATA  ACQUISITION 

The  LOGACS  data  were  tape  recorded  on-orbit  and  then,  at  the 
ground  station  acquisitions,  the  tape  recorder  was  commanded  to  reproduce. 
The  piayback  of  the  orbital  data  was  recorded  at  the  ground  station.  It  should 
be  made  clear  that  while  the  spacecraft  tape  recorder  was  in  reproduce  mode, 
the  LOGACS  data  were  also  being  transmitted  in  real  time  on  another  telem¬ 
etry  channel.  Thus,  no  information  was  lost  during  the  tape  recorder 
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reproduce  sequence.  The  real-time  data  were  decommutated  and  printed  out 
subsequent  to  each  station  acquisition.  Analog  records  of  much  of  the  real¬ 
time  data  were  also  made  during  many  station  acquisitions. 

After  the  LOGACS  flight,  all  the  recordings  mace  at  the  ground 
stations  were  collected.  These  included  both  the  real-time  data  and  the  tape 
recorder  playback  data.  All  of  this  information  was  decommutated,  changed 
to  desired  format,  ordered  with  respect  to  the  experiment  clock  time,  and 
encoded  on  a  floating  point  binary  tape  suitable  for  computer  processing. 
Listings  of  the  final  merged  tape  were  made,  as  were  plots  of  the  MESA  out¬ 
put  pulse  rate.  The  attitude  control  jet  firing  occurrences  were  plotted  on 
the  same  graph  as  the  MESA  output  pulse  rate,  as  were  the  Agena  attitude 
control  system  gyro  signals.  A  sample  of  these  plots  and  typical  data  pages 
are  shown  in  Fig.  F-l  and  Table  2. 

The  second  channel  of  the  tape  recorder  contained  the  formated 
accelerometer  output  pulses.  The  information  from  this  channel  has  not  been 
exploited  because  the  accelerometer  output  pulses  per  second  were  available 
on  the  pair  of  counters  that  were  on  the  commutated  channel.  While  the  pulse 
information  on  the  second  channel  generally  duplicates  the  counter  informa¬ 
tion,  it  does  contain  some  fine  grain  information  which  is  not  apparent  from 
the  accelerometer  output  pulse  counters. 

I.  DATA  PROCESSING 

The  objective  of  the  LOGACS  data  processing  was  to  determine 
values  of  MESA  instrument  bias,  instrument  scale  factor,  and  drag  on  the 
vehicle.  The  specifications  for  the  LOGACS  data  analysis  program  are 
delineated  in  Appendix  C.  A  listing  of  the  LOGACS  de  .a  analysis  program 
is  presented  in  Appendix  D. 
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The  essence  of  the  LOGACS  data  processing  was  to  fit,  in  a  least 
squares  sense,  the  data  to  the  following  simple  model 


A  =  K 
m 


B  I  Rw  + 


VAdrag  +  Adragf  +  Adrag  T  COS  (wt  +  *> 


where 

A^  =  accelerometer  output  (pulses /sec) 

B  =  instrument  bias  (g) 

K  =  instrument  scale  factor  (pulses / sec -g) 

R  =  distance  from  centrifuge  center  of  rotation  to  center  of  mass 
of  the  MESA  proof  mass 

u  =  centrifuge  rotation  speed  (rad/sec) 

4>  =  phase  angle  of  centrifuge  (rad) 

A,  =  vehicle  acceleration  due  to  noncons  ergative  forces  acting  on 
C  ra^  the  vehicle  (g) 

The  values  of  B,  K,  A^  ,  and  <f>  are  solved  for  during  each  mode  sequence; 
that  is,  a  solution  for  the  above  parameters  is  accomplished  for  every  2048- 
sec  period  during  which  no  large  data  gap  occurs. 

An  alternate  data  processing  routine  was  also  employed  to  independ¬ 
ently  solve  for  the  same  parameters.  The  alternate  scheme  was  conceptually 
almost  identical  to  the  above  scheme  but  did  differ  in  detail.  A  comprehensive 
description  of  the  alternate  processing  scheme  is  contained  in  Section  IL 

J.  DISCUSSION  OF  RESULTS 

1.  INSTRUMENT  RESULTS 

The  MESA  instrument  data  obtained  from  the  data  analysis  program 
that  made  use  of  a  second -order  fit  to  the  LOGACS  Mode  1  and  Mode  2  data 
(see  Appendix  Dl  are  presented  in  Fables  3  and  4  and  Figs.  13  and  14. 

Table  3  shows  the  instrument  bias  and  scale  factor  where  no  edit  was  made 
of  thrust  valve  firings;  Fig.  13  is  a  plot  of  the  data  in  Fable  3.  Table  4  shows 
the  bias  and  scale  factor  for  which  an  edit  was  made  of  all  data  periods  con¬ 
taining  thrust  valve  firings;  Fig.  14  is  a  plot  of  the  data  in  Fable  4. 
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La ii-M  i.  MES/*  Instrument  Bias  and  Scale  r at  tor 
(no  edit  uf  data  acquired  during  thrust 
valve  firings! 


Experiment 

Cli  ck  rim-3, 

SeC 

Bias, ^ 

R  X  1C  ' 

j 

! 

Scale  F actor , 
pulses/ sec -g  X  30' 

12  288 

i  -  >.  s 

0. 9860 

id  6  62  4 

i  -  3.  0 

0. 9863 

28672 

-  3  Q 

*  ' 

0.  983.3 

30720 

i  1  °*  5 

0.  9840 

4?  056 

-  t>.  6 

0. 9856 

47104 

-  2.8 

0.  9891 

4  H2 

-  5.  9 

0. 9848 

-•1200 

6.  7 

0. 9857 

"•  1248 

0. 9815 

5  7296 

-  7) .  8 

0. 9853 

5  7844 

8.  8 

0. 9919 

5  0  3"  2 

-  7.  < 

0. 9845 

t  3488 

*  2. 

0. 9795 

'.55  <8 

6.  2 

0. 9853 

i  '8  84 

-  4.3 

0. 9825 

6  u32 

r  .  3 

0. 9836 

71680 

'  .  •' 

i  0.9853 

1  n  ,'4 )  ( 

—  .  , 

0. 9843 

110404 

r. .  i 

0. 9852 

11  2">  12 

0.  9843 

1 14460 

i  8.  t. 

0. 9852 

126848 

1  0.  8 

0. 9891 

128896 

i  *  3.  4 

0.9836 

130944 

!  ■  7.  5 

0. 9835 

!  12  >92 

-  =.  1 

0. 9838 

195520 

-  7 .  7 

0. 9833 

1575.0 

-  7.  9 

0.  9830 

i  6 1664 

■  7.  9 

0. 9834 

103712 

,  } 

0. 9824 

173824 

,  2 

0.  '"'851 

175872 

4.  7 

0. 9853 

1  7  7920 

-  8.  5 

0. 9846 

179968 

i  '» 

-  < .  - 

0. 9825 

186112 

1  0.  9 

0, 9972 

188160 

-  7.  8 

0.  ‘>842 

194304 

-  8.4 

0. 984  3 

1  "53 


Tabic  3.  MESA  Instrument  Bias  and  Scale  Fact'" 
(no  edit  of  data  acquired  during  thrust 
valve  firings)  (continued) 


Experiment 

Clock  Time, 
sec 

Bias, 
g  X  10-7 

Scale  F actor, 
pulse s/ sec  -g  X  1 07 

196352 

-  6.3 

0.9799 

198400 

-  9.4 

0.  9859 

204544 

-  6.9 

0. 9827 

208640 

-  0.2 

0.9776 

210688 

-  8.  1 

0. 9841 

214784 

-  7.7 

0. 9837 

218880 

-  9.  1 

0. 9904 

220928 

-  8.  3 

0. 9825 

229120 

-12.6 

0.9899 

231168 

-  9.4 

0. 9850 

241408 

-  8.  3 

0.9821 

243456 

-  8.8 

0. 9837 

247552 

-  9.2 

0. 9842 

249600 

-  8.1 

0.9758  i 

251648 

0.  0 

0. 9743 

253696 

-  8.  5 

0. 9821 

257792 

+  73.  4 

0. 9849 

259712 

+  72.  8 

0.9848 

269952 

)  7  3.  0 

0.9847 

272000 

+  75,  6 

0.9913 

274048 

+  73.  9 

0.  9835 

280192 

+  72.  2 

C. 9846 

284288 

+  71.7 

0.9868 

286336 

+  71.  6 

0.9846 

288384 

+66.  8 

0. 9940 

290434 

+  72.6 

0.  9838 

204528 

-72.  5 

0.9842 

296576 

+  71.  3 

0. 9866 

298624 

-12.  5 

1.0193 

300672 

-  6.4 

0. 9840 

302720 

-  6.9 

0.9830 

304768 

+  21.  B 

0. 9563 

306816 

-  6.  8 

0.9828 

317056 

-  7.4 

0. 9843 

337536 

-  6.  1 

0.9815 

j  339584 

-  8.  3 

0. 9835 

349696 

-  8.  6 

0. 9836 

355840 

-  8.5 

0.  9821 

NOTE:  The  value  of  1  g  =  32.  145  ft/ see 2  was  used  to  calculate  the  above  data 
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Table  4.  MESA  Instrument  Bias  and  Scale  Factor 
(edit  of  all  data  acquired  during  thrust 
valve  firings)  (continued) 


Experiment 

Clock  Time,” 
sec 

Bias, 
g  X  10  -7 

Scale  F actor, 
pulses/sec-g  X  10^ 

247552 

-  9.2 

0.9835 

249600 

-  7.  1 

0.  9759 

253696 

-  9.2 

0.9824 

257792 

+  73.  0 

0.9851 

259712 

+  74.  5 

0.9827 

269952 

+  72.  5 

0.9850 

274048 

+  71.7 

0.9851 

280192 

+  71.  8 

0. 9848 

284288 

+  67.4 

0.9908 

286336 

+  71.9 

0.9843 

i  290434 

+  71.4 

0. 9848 

296576 

+  70.  9 

0.9867 

1  300672 

-  8,2 

0.9858 

302720 

-  6.6 

0.9709 

306816 

-  7.0 

0.9831 

317056 

-  8.7 

0.9860 

337536 

-15.7 

0.9907 

339584 

-  8.4 

0.9834 

349696 

-  8.7 

0.9835 

355840 

-10.  7 

0.9834 

The  time  sequence  in  this  table  differs  from  that  of  Table  3;  points  missing 
here  were  edited  because  of  thrust  valve  firings. 

NOTE:  The  value  of  1  g  ~  32.  145  ft/sec^  was  used  to  calculate  the  above 
data. 

The  MESA  instrument  data  obtained  from  the  alternate  data  analysis 
technique  differed  from  the  second-order  fit  data  analysis  in  that  the  editing 
was  done  by  hand  and,  therefore,  involved  judgment.  Table  5  shows  the  bias 
and  scale  factor  values  determined  by  this  alternate  method,  and  Fig,  15  is  a 
plot  of  the  data  in  Table  5. 

A  summary  of  the  mean  values  and  the  standard  deviation  about  the 
mean  for  the  MESA  instrument  bias  and  scale  factor  is  presented  in  Table  6. 
The  values  of  bias  and  scale  factor  presented  in  Table  6  are  derived  from  that 
set  of  data  points  which  was  common  to  the  three  different  analysis  programs. 
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Data  Analysis  with  No  Thruster  Edit 
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ig.  14.  Data  Analysis  with  Thruster  Edit 


Table  5.  MESA  Instrument  Bias  and  Scale  Factor 
(data  reduced  by  means  of  alternate 
analysis  program) 


Experiment 

Clock  Time,5' 
sec 

Bias, 

g  xio-7 

Scale  Factor, 
pulses/sec-g  X  107 

12288 

-  4.  1 

0.9861 

28672 

-  4.2 

0.9834 

32768 

-  5.0 

0.9840 

45056 

-  6.5 

0.9855 

49152 

-  6.6 

0.9852 

51200 

-  6.4 

0.9843 

55296 

-  6.9 

0.9861 

59392 

-  6.9 

0.9853 

65536 

-  7.2 

0.9853 

71680 

-  7.0 

0.9844 

108416 

-  7.8 

0.9841 

114560 

-  8.0 

0.9849 

124800 

-  8.9 

0.9840 

128296 

-  8.6 

0.9838 

130944 

-  8.4 

0.9848 

135040 

-  8.1 

0.9832 

157568 

-  8.3 

0.9838 

161664 

-  8.0 

0.9832 

163712 

-  7.5 

0.9817 

173824 

-  8.8 

0.9847 

177920 

-  8.4 

0.9845 

179968 

-  7.8 

0.9828 

188160 

-  8.4 

0. 9844 

190208 

-  7.7 

0.9842 

194304 

-  8.  0 

0.9833 

196352 

-  8.0 

0.9815 

200448 

-  8.6 

0.9850 

204544 

-  7.9 

0.9834 

210688 

-  8.6 

0.9846 

214784 

-  9.1 

0.9852 

216832 

-  7.4 

0.9824 

220928 

-  8.3 

0.9823 

227072 

-  8.8 

0.9837 

231168 

-  9.7 

0.9847 

233216 

-  8.7 

0.9828 

243456 

-  9.0 

0.9837 

The  time  sequence  in  this  table  is  different  from  that  of  Table  3;  points 
missing  here  were  edited  because  of  thrust  valve  firings. 


Tabie  5.  MESA  Instrument  Bias  and  Scale  Factor 
(data  reduced  by  means  of  alternate 
analysis  program)  (continued) 


Experiment 

Clock  Time,* 
sec 

Bias, 
g  X  10-7 

Scale  F actor, 
pulses/sec-g  X  10? 

247552 

-  9.2 

0.9838 

249600 

-  8.0 

0.9802 

253696 

-  8.8 

0.9829 

257792 

+73.4 

0.9847 

259712 

+74.  1 

0.9836 

274048 

+72.9 

0.9846 

276096 

+71.9 

0.9848 

280192 

+72.2 

0. 9844 

286336 

+72.  2 

0.9839 

290432 

+72.  1 

0.9841 

296576 

+70.  8 

0.9869 

300672 

-  7.5 

0.9848 

302720 

-  7.6 

0.9842 

306816 

-  7.4 

0.9835 

312960 

-  8.0 

0.9850 

317056 

-  8.6 

0.9856 

329344 

-  8.0 

0.9830 

333440 

-  8.1 

0.9835 

339584 

-  8.5 

0.9841 

343680 

-  8.5 

0.9835 

349696 

-  8.3 

0.9831 

355840 

-  8.7 

0.9821 

The  time  aeouence  in  this  table  differs  from  that  of  Table  3;  points  missing 
here  were  edited  because  of  thrust  valve  firings. 

A  ? 

NOTE:  The  value  of  1  g  —  32.  145  ft/sec  was  used  to  calculate  the  above 
data. 


The  decision  to  employ  only  the  common  data  points  was  made  in  an  effort  to 
avoid  the  errors  introduced  by  the  computational  procedures.  While  the  raw 
input  data  to  the  three  analysis  programs  were  identical,  the  degree  and 
specific  manner  in  which  the  edit  of  the  data  was  performed  did  influence 
some  of  the  bias  and  scale  factor  determinations  in  unique  ways. 

The  stability  of  the  bias  and  scale  factor  during  the  nearly  100-hr 
orbital  calibration  is  evident  from  the  small  values  of  standard  deviation 
shown  in  Table  6  as  well  as  from  the  data  plots  shown  in  Figs.  13  through  15. 
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iate  Data  Analysis 


The  relationship  between  the  values  of  bias  and  scale  factor 
determined  prior  to  flight,  on-orbit  with  the  7.  5  X  10*  g  cross -axis  support, 
and  on-orbit  with  the  1.  5  g  cross-axis  support  are  presented  in  Table  7. 


Table  7.  Expected  Relationship  of  Bias  and  Scale  Factor 


Preflight 

On-Orbit 

7.  5  X  10"^  g  Support 

On-Orbit 

1. 5  g  Support 

Scale  factor 

Bias 

X 

Y 

V 

-'V 

(1/200)  Y 

X 

Y 

The  bias  and  scale  factor  coefficients  for  the  MESA  instrument  were 
determined  by  Bell  Aerosystems  Company  prior  to  installation  of  the  MESA 
in  the  centrifuge  assembly.  The  uncertainty  in  the  value  of  bias  as  measured 

_7 

at  Bell  was  on  the  order  ±50  x  10  g.  The  uncertainty  in  the  value  of  scale 
factor  as  measured  at  Bell  was  on  the  order  ±0.5  percent. 

The  bias  was  measured  at  Bell  on  three  separate  occasions.  The 
result  of  the  three  determinations  is  presented  in  Table  8. 

Table  8.  Bias  Determination  at  Bell  Aerosystems  Company 

.  7 

Determination  No.  1  0  ±  50  X  10  g 

.7 

Determination  No.  2  0  ±  50  X  10  g 

Determination  No.  3  +50  ±  50  X  10"^  g 

The  scale  factor  was  measured  once  at  Bell  and  the  value  was  found 

7 

to  be  1.  031  ±  0.  005  X  10  pulses/sec  -g. 

If  the  values  of  bias  and  scale  factor  obtained  prior  to  the  flight  are 
compared  to  those  obtained  on-orbit,  we  find  two  startling  facts.  The  first 
is  that  the  scale  factor  changed  by  about  4.  7  percent,  although  no  change  was 
expected.  The  second  is  that  the  preflight  bias  value  correlates  reasonably 
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well  with  the  on-orbit  v:.lue  when  the  1.  5  g  cross -axis  support  was  used,  but 
when  the  cross -axis  support  voltages  were  reduced  by  a  factor  of  200,  the 
bias  changed  polarity  from  positive  to  negative  and  became  a  factor  of  10 
smaller. 

Since  the  LOGACS  flight,  the  above  anomalies  in  bias  and  scale 
factor  have  been  given  a  considerable  amount  of  attention.  In  the  case  of  the 
bias  changing  in  a  manner  not  linearly  correlated  with  the  cross -axis  support 
voltage,  a  clear  explanation  has  evolved  as  to  why  the  bias  should  in  general 
never  be  expected  to  correlate  linearly  with  the  cross -axis  support  voltages. 

Consider  the  cross-axis  support  forces  that  act  upon  the  proof  mass 
of  the  accelerometer.  For  ease  of  illustration  let  us  confine  our  examination 
to  a  single  plane.  The  forces  on  the  proof  mass  when  the  proof  mass  does 
not  possess  symmetry  would  be  as  in  Fig.  16. 


Fig.  16.  Nonsymmetrical  Proof  Mass 

The  individual  forces  on  the  proof  mass  are  the  results  of  the  electro¬ 
static  force  assoc. ated  with  each  of  the  eight  support  electrodes.  Each  support 
force  acts  normal  to  the  region  of  the  proof  mass  which  it  is  intended  to  sup¬ 
port.  Thus,  the  forces  F^  and  F^  are  acting  normal  to  the  sensitive  axis 
direction  for  the  geometry  shown  in  Fig.  16,  while  forces  F^  and  are 
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acting  at  an  angle  a  away  from  such  a  normal.  Clearly  then,  the  force  on 
the  proof  mass  in  the  direction  of  the  sensitive  axis  is  the  instrument  bias. 
The  bias  can  be  expressed  as 

Bias  =  F3sina  +  F4sinc'  (10) 

or 

Bias  =  (F ^  +  F^)  sin  a  (11) 

Also,  the  support  forces  must  satisfy  the  relationship 

Fj  +  Fz  =  (F3  +  F4)  cos  a  (12) 

if  static  equilibrium  is  to  be  maintained  in  the  vertical  direction.  So  far,  we 
have  not  taken  into  account  the  proof  mass  itself  and  the  alterations  in  the 
support  forces  that  occur  as  a  function  of  the  acceleration  in  the  cross -axis 
direction;  that  is,  the  relationships  above  are  only  valid  for  a  zero  cross - 
axis  acceleration  condition. 

Now,  imagine  that  the  conditions  imposed  above  for  Fig.  16  are 
changed  such  that  a  cross -axis  acceleration,  a,  exists  and  is  directed 
vertically  upward.  The  imposition  of  the  cross-axis  acceleration  changes 
the  equilibrium  in  the  vertical  direction  by  the  force  ma.  Thus,  we  have 

Fj+F2=Fj+F2*  ma  (13) 

and 

(F^  +  F 4>  cos  a  =  (Fj  +  F^)  cos  a  +  ma  (14) 
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where  the  prime  notation  signifies  the  condition  existing  with  the  cross -axis 
acceleration.  For  this  condition,  we  find  that 

Bias  =  (Fj  +  F^)  sin  a  (15) 

Substituting  Eq.  (14)  into  Eq,  (15),  we  have 

Bias  =  (Fj  +  F^)  sin  a  +  ma  tan  a  (16) 

Thus,  the  bias  has  changed  from  the  value  shown  in  Eq.  (11)  by  the  simple 
expedient  of  changing  the  cross -axis  acceleration.  The  same  general  argu¬ 
ment  can  be  extended  to  show  that  any  asymmetry  in  the  proof  mass  or  sup¬ 
port  electrode  geometry  will  lead  to  a  cross -axis  acceleration  sensitivity  of 
the  bias. 

The  net  cross -axis  acceleration  sensitivity  for  the  MESA  used  in 
LOGACS  was  on  the  order  of  50  pg/cross-axis  g.  Recent  analysis  by  the 
MESA  manufacturer  (Bell  Aerosystems  Company!  suggests  that  the  sensitivity 
can  be  reduced  to  the  order  of  5  pg/cross-axis  g. 

The  reason  for  the  scale  factor  change  of  4.7  percent  is  not  specifi¬ 
cally  known.  The  source  of  a  change  in  scale  factor  can  be  generally  sepa¬ 
rated  into  electrical  or  mechanical  categories. 

The  mechanical  category  includes:  (a)  the  length  of  the  radius  arm 
of  the  centrifuge  being  different  from  that  used  to  compute  the  orbital  scale 
factor  values,  (b)  the  angular  velocity  of  the  centrifuge  being  different  from 
that  used  to  compute  the  orbital  scale  factor  values,  and  (c)  the  internal 
geometry  of  the  MESA  shifting  after  the  ground  calibration. 

The  radius  arm  length  computations  were  carefully  rechecked  after 
the  flight  to  ascertain  whether  an  error  could  have  been  made.  Since  the 
nominal  length  of  the  radius  arms  is  4  in.,  the  necessary  change  in  radius 
arm  length  to  produce  a  4.7  percent  change  in  scale  factor  would  have  to  be 
about  0.  2  in.  Such  a  gross  error  could  only  result  from  something  like  a 
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transposition  of  numbers  or  a  computational  oversight.  The  rechecking  of 
the  radius  arm  computations  failed  to  uncover  any  such  error;  hence  we  can 
discard  a  change  in  radius  arm  length  as  the  source  of  the  scale  factor  shift. 

In  the  course  of  the  meticulous  rechecking  of  the  radius  arm  length, 

no  reason  was  found  to  douot  the  measured  value  of  R  of  3.  9701  ±  0.0024  in. 

The  uncertainty  in  centrifuge  rotational  speed  (w)  was  less  than  20  parts  per 

million  as  judged  by  the  demonstrated  accuracy  of  the  experiment  clock,  and 

2 

therefore  it  can  be  neglected  insofar  as  its  contribution  to  uncertainty  in  Rw 
is  concerned.  Thus,  the  uncertainty  in  scale  factor  due  to  uncertainties  in  R 
and  ijj  is  due  almost  completely  to  the  uncertainty  in  R,  which  leads  to  ±0.06 
percent  uncertainty  in  the  absolute  scale  factor.  This  means  that  the  mean 
value  of  scale  factor  determined  from  the  flight  data  is  uncertain  in  absolute 
value  by  ±0.06  percent. 

The  rotational  speed  of  the  centrifuge  was  independently  checked  by 
means  of  the  telemetry  data.  Each  rev.  of  the  rotating  arms  of  the 
centrifuge  tripped  a  pair  of  reed  relay  position  pickoffs  located  180  deg  apart. 
Thus,  obtaining  the  time  between  the  tripping  of  the  reed  relays  from  the 
real  time  telemetry  data  made  it  possible  to  show  that  the  rotational  speed 
was  at  the  design  speed  with  an  error  not  exceeding  0.  1  percent  and  not  at 
the  speed  required  to  produce  a  4.  7  percent  scale  factor  error. 

The  possibility  of  a  mechanical  shift  of  elements  within  the  MESA 
exists  but  is  considered  to  be  an  unlikely  source.  The  demonstrated  stability 
of  scale  factor  during  the  100 -hr  flight  serves  to  testify  that  a  gradual  shift 
of  the  mechanical  elements  of  the  MESA  did  not  occur.  The  remaining  possi¬ 
bility  is  that  a  shift  of  the  mechanical  elements  occurred  as  a  step  change  and 
that  thereafter  the  elements  remained  fixed.  The  latter  does  exist  as  a  finite 
possibility,  but  it  appears  contrived  and  is  not  considered  to  be  a  likely 
explanation  for  the  4.7  peveent  scale  factor  shift. 

Another  possible  explanation  for  the  scale  factor  shift  is  that  it  does 
not  exist  at  all  but  is  simply  the  result  of  an  erroneous  ground  determination. 
This  possibility  fails  because  independent  sources  have  examined  the  data 
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and  calculated  the  scale  factor  only  to  verify  that  a  4.  7  percent  shift  in  scale 
factor  actually  does  exist. 

The  possible  sources  of  scale  factor  shift  which  fall  into  the  electri¬ 
cal  category  involve  two  basic  circuits:  (a)  the  temperature  controller  for 
the  magnetic  core  that  supplies  the  precision  rebalance  pulses  and  (b)  the 
capacitance  pickcff  Arhich  detects  the  proof  mass  displacement  along  the 
sensitive  axis  and  which  is  used  to  trigger  the  rebalance  pulses. 

The  change  in  set  point  or  control  temperature  which  would  produce 
a  change  in  pulse  size  (i.  e.,  scale  factor)  of  4.7  percent  is  approximately 
94°F.  Such  a  change  is  considered  extremely  unlikely  because  the  orbital 
data  are  stable  during  the  entire  100-hr  flight.  Again  the  possibility  exists 
that  a  step  change  in  set  point  of  the  temperature  controller  occurred  after 
the  ground  test  but  remained  stable  afterwards.  This  latter  possibility  again 
appears  contrived  and  is  not  considered  to  be  a  likely  explanation  for  the 
4.  7  percent  scale  factor  shift. 

A  change  in  scale  factor  will  result  from  a  shift  in  the  electrical  null 
of  the  capacitance  bridge  pickoff.  Figure  17  illustrates  how  the  scale  factor 
changes  as  a  function  of  shifts  in  the  electrical  null  of  the  pickoff.  The 
figure  is  not  to  scale  but  is  used  to  illustrate  the  principles  involved,  As  the 
proof  mass  moves  along  the  sensitive  axis,  the  capacitance  bridge  becomes 
unbalanced  and  the  amplified  and  phase -demodulated  output  of  the  bridge  is 
directly  proportional  to  its  displacement.  The  bridge  output  voltage  is  repre¬ 
sented  as  E  in  Fig.  17.  Imagine  that  the  pickoff  electrical  null  is  represented 
by  Eq  and  that  E^,  and  -E^,  represent,  respectively,  the  positive  and  negative 
thresholds  that  are  used  to  gate  the  voltage  pulses  to  the  forcer  rings  and 
thus  apply  force  to  the  proof  mass.  Assume  that  the  scale  factor  for  each 
sense  is  identical  for  this  case.  The  force  on  the  proof  mass  during  each 
application  of  voltage  to  the  forcer  rings  is 

force  in  positive  sense  oc  l /A*" 

force  in  negative  sense  oc  l/A*" 
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I-' iff.  17.  Relationship  of  Scale  Factor  to  Pickoff  Null 

Now  imagine  that  the  electrical  null  of  the  bridge  shifts  (i.  e.  ,  one  leg  of  the 
bridge  changes  in  capacitance).  ’1  hen  the  proof  mass  must  move  to  a  new 
position  in  order  for  the  bridge  to  be  balanced.  The  new  null  position  of  the 
proof  mass  is  represented  in  Fig.  17  as  EfJ.  The  proof  mass  displacement 
corresponding  to  the  shift  in  electrical  null  is  represented  as  6 ,  Due  to  the 
shift  in  balance  of  the  capacitance  bridge,  we  find  that  the  proof  mass  posi¬ 
tions  corresponding  to  the  positive  and  negative  thresholds  have  changed.  The 

f  t 

now  threshold  positions  are  represented  by  E.^  and  It  should  be  noted 

that  E.,,  and  -E.^  have  the  same  magnitude  as  E.^  and  -  E^.  The  prime  nota- 
lion  is  used  as  an  expedient  to  denote  the  change  in  position  which  is  due  to 
the  change  in  balance  of  the  bridge.  In  this  case,  the  force  on  the  proof  mass 
during  rath  application  of  voltage  to  the  forcer  rings  is 


(force)  in  positive  s-  nsc  <r. 


(> « f") 


1  1 

(force)  m  negative  sense  no - y  — p — - - 

f  A  t  >■,  £  <■  l<  A 
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Since  the  force  applied  to  the  proof  mass  during  the  application  of  a  voltage 
pulse  to  '-he  forcer  rings  is  different  when  the  bridge  balance  is  changed,  the 
pulse  rate  that  must  be  supplied  to  the  forcer  rings  to  rebalance  the  proof 
mass  for  a  given  acceleration  will  change  in  the  same  manner. 

The  capacitance  bridge  unbalance  necessary  to  produce  a  scale 
factor  change  of  4.7  percent  is  very  small  --on  the  order  of  0.070  pF. 

(Since  the  pickoff  gain  is  on  the  order  of  1000  pF/cm,  the  displacement  6  that 

_5 

corresponds  to  the  0.  070-pF  bridge  unbalance  is  about  7  X  10  cm.) 

While  the  precise  mechanism  for  a  shift  in  the  capacitance  bridge 
balance  is  not  known,  it  appears  to  be  the  most  likely  explanation  for  a  large 
change  in  scale  factor. 

The  expressions  developed  above  for  the  force  applied  to  the  proof 
mass  when  the  pickoff  null  has  been  displaced  reveal  that  a  scale  factor  shift 
caused  in  this  manner  would  exhibit  a  differential  behavior;  that  is,  the  scale 
factor  in  one  sense  would  increase  by  1  +  26/A  whereas  the  scale  factor  for 
the  opposite  sense  would  decrease  by  1  -  26/A.  Recognition  of  this  led  to 
the  hope  that  such  a  behavior  could  be  recognized  in  the  L.OGACS  flight  data 
and  hence  lend  support  to  the  conjecture  that  a  shift  in  the  capacitance  bridge 
balance  caused  the  scale  factor  shift. 

The  portion  of  the  data  that  lends  itself  to  testing  for  a  different 
scale  factor  in  the  two  senses  is  where  mode  3  (MESA  pointing  forward)  ends. 
The  centrifuge  then  rotates  180  deg  and  mode  4  (MESA  pointing  aft)  begins. 
Two  things  interfere  with  any  attempt  to  examine  the  data  for  differences  in 
scale  factor  which  are  associated  with  sense.  Clearly  a  large  drag  value  is 
desirable  to  test  for  the  scale  factor  change,  but  it  is  also  true  that  the  drag 
value  changes  more  rapidly  as  the  drag  value  gets  large.  Another  associated 
difficulty  is  that  the  vehicle  attitude  control  firings  become  significantly  more 
,'requent  as  drag  increases  and  hence  the  "noise"  in  the  data  induced  by  the 
attitude  control  firings  increases.  Thus,  we  find  ourselves  faced  with  the 
unfortunate  conflict  wherein  the  drag  and  the  "noise"  on  the  drag  signal 
increase  substantially  together.  As  a  result,  the  attempts  to  date  to  show 
that  a  differential  relationship  exists  for  the  scale  factor  for  each  sense  have 
failed. 
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The  MESA  manufacturer  also  produces  a  line  of  accelerometers 
that  employ  a  hinged  pendulum  as  the  proof  mass.  These  accelerometers 
make  use  of  a  very  sensitive  capacitance  bridge  pickoff  (=600  pF/cm)  which 
is  similar  in  both  concept  and  design  to  that  of  the  MESA.  Due  to  the 
mechanization  of  the  pendulous  accelerometers,  a  shift  in  the  balance  of  the 
capacitance  bridge  pickoff  results  in  a  shift  in  the  null  bias  of  the  instrument. 

A  change  of  0.  07  pF  in  the  balance  of  the  capacitance  bridge  of  the  pendulous 
accelerometer  would  result  in  a  null  bias  change  of  350  ng.  Such  a  shift 
would  be  abnormally  large  for  these  devices.  The  similarity  in  the  capaci¬ 
tance  pickoffs  would  suggest  that  if  a  shift  of  0.07  pF  is  unusual  in  the  hinged 
pendulum  accelerometers,  it  would  also  be  unusual  in  the  MESA.  Further¬ 
more,  since  the  MESA  scale  factor  demonstrated  approximately  0.  1  percent 
stability  during  100  hr  on-orbit,  the  shift  would  have  had  to  be  a  step  change 
rather  than  a  random  wander. 

The  facts  above  do  not  prove  but  do  lend  support  to  the  following  con¬ 
jecture:  if  the  scale  factor  shift  observed  in  the  LOGACS  MESA  was  due  to  a 
shift  in  the  balance  of  the  capacitance  bridge,  then  that  shift  was  an  anomalous 
one  and  should  not  be  considered  typical  of  the  MESA  instrument. 

In  Section  G  it  was  noted  that  the  accelerometer  remained  saturated 
for  a  long  period  of  time  after  launch  and  again  after  the  orbit  adjust  was 
accomplished.  The  reason  for  the  long  saturation  time  has  not  been  con¬ 
clusively  proved,  but  a  plausible  explanation  has  evolved. 

When  the  MESA  proof  mass  is  driven  over  against  one  side  of  the 

Instrument,  the  force  which  can  be  applied  to  restore  the  proof  mass  to  the 

centered  position  is  only  1/4  that  which  can  be  applied  when  the  instrument 

-4 

is  centered.  Thus,  with  a  1-gram  proof  mass  and  a  5  x  10  g  constrainment 
capability  along  the  sensitive  axis,  the  force  that  can  be  applied  to  the  centered 
proof  mass  is  about  5  X10  dyne.  If  the  proof  mass  moves  against  the 
inactive  forcer  ring  set  (and  hence  is  twice  as  far  away  from  the  operating 
forcer  rings),  the  force  that  can  be  applied  is  about  1.25  X  10  dyne.  Now, 
if  the  average  force  on  the  proof  mass  approaches  1.  25  X  10  *  dyne  and  if  the 
proof  mass  is  already  over  against  the  stop,  then  the  proof  mass  will  take 


longer  and  longer  to  return  to  the  centered  position.  In  fact,  il  the  force 
that  can  be  applied  exactly  equals  or  is  less  than  the  average  force  on  the 
proof  mass,  then  the  proof  mass  will  remain  against  the  stop  indefinitely. 
Laboratory  tests  with  a  MESA  instrument  have  demonstrated  that  long 
saturation  times  can  be  easily  produced  by  the  technique  described  above. 
A  typical  curve  generated  in  the  course  of  such  tests  is  shown  in  Fig.  18. 
While  the  tests  do  not  show  conclusively  that  this  phenomenon  explains  the 
observed  behavior,  they  lend  strong  support  to  this  hypothesis. 
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TEST  procedure: 

PROOF  MASS  IS  PERMITTED  TO  TRAVEL 
OVER  TO  A  STOP.  SENSITIVE  AXIS  IS 
THEN  ORIENTED  TO  SOME  POSITION  WHICH 
INDUCES  AN  ACCELERATION  EQUIVALENT 
TO  A  FIXEO  PERCENT  OF  FULL  SCALE. 
RECOVERY  TIME  IS  CALCULATED  FROM 
THE  TIME  THE  SENSITIVE  AXIS  IS 
ORIENTED  UNTIL  THE  ACCELEROMETER 
OUTPUT  RETURNS  TO  AN  ON  -SCALE 
READING 


5  10  15  20  25 

ACCELERATION  APPLIED  ALONG  SENSITIVE  AXIS,  PERCENT  OF  FULL  SCALE 


Fig.  18.  Typical  Curve  from  MESA  Saturation  Tests 
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SECTION  II 


LOGACS  DATA  REDUCTION 

A.  LOGACS  SENSING  ENVIRONMENT 

Since  it  is  desired  to  calibrate  the  accelerometer  in  a  drag-free 
environment  while  it  is  desired  to  measure  drag  with  the  accelerometer  in  a 
fixed  orientation,  the  environment  for  either  objective  was  not  ideal.  In  addi¬ 
tion,  during  rotating  modes  the  lateral  forces  and  moments  induced  by  aero¬ 
dynamics  and  the  Agena  control  system  contribute  to  the  measured  accelera¬ 
tion.  For  these  reasons  it  was  necessary  to  perform  an  analysis  of  the 
sensing  environment.  This  section  presents  this  analysis  and  develops  the 
algorithm  lor  the  acceleration  sensed  by  the  LOGACS  accelerometer.  The 
terms  in  this  equation  are  then  analyzed  and  their  order  of  magnitude  estab¬ 
lished.  Based  on  this  analysis,  simplifications  are  made  to  the  algorithm 
used  for  processing  the  accelerometer  data. 

1.  AGENA  CONFIGURATION 

Figure  19  presents  a  schematic  configuration  diagram  of  the  vehicle 
and  the  LOGACS  instrument.  A  body-fixed  coordinate  system  is  used  as  a 
reference  for  defining  the  terms  to  be  used  in  the  analysis  which  follows. 

Th  x  (roll)  axis  lies  on  the  geometric  center  line  of  the  vehicle  pointing 
forward,  the  y  (pitch)  axis  points  out  the  right  side  of  vehicle,  and  the 
z  (yaw)  axis  is  down.  The  center  of  this  coordinate  system  is  defined  to  lie 
on  the  vehicle  center  line  at  the  reported  Agena  eg  station  (251.4  in.).  Based 
on  these  definitions  and  those  in  Fig.  19,  the  data  pertinent  to  the  configura¬ 
tion  used  in  this  analysis  are  presented  in  Table  9. 
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Table  9.  Configuration  Data 


Term 

Symbol 

Value 

Units 

Rema  rks 

Mass 

m 

59.  5 

slug 

1912  !b 

‘x 

1™ 

slug-ft^ 

Inertia 

1 

V 

2560 

15  lug  -  ft^ 

! 

z 

2560 

slug  -  ft^ 

X 

<•« 

0 

ft 

Station  351.4  in. 

Center  of  gravity 

0.  1*15 

ft 

2,  3  in. 

Zog 

0.  125 

ft 

1.  5  in. 

X 

ac 

.<» 

ft 

i.OGACS  location 

Y 

*  ac 

(J 

ft 

Z 

ac 

X 

1.  5 

.  '< 

ft 

Gas  valve  locations 

' 

Y 

«v 

'1 

ft 

z 

gv 

*2.  If  6 

ft 

X 

CP 

See  Section  2.  5 

Aerodynamic  center 

Y 

CP 

0 

ft 

i 

/. 

cp 

0 

11 

Reference  area 

s 

1 ".  6 

fta 

Reference  length 

'-r 

5 

ft 

]  Gas  valve  force 

r. 

i 

0.  5 

lb 

(i  -  1 1  2.t  •  •  •  t  6l 

j 

;  LCjOACS  radius  arm 

r 

<t 

1/9 

ft 

|  Turntable  phase  angle 

- 180  s  ♦  S  180 

deg 

} 

1 

0 

r  pm 

Mode  3(«t  --  0),  Mode  4(4(  »  *> 

!  Turntable  rates 

6.445 

rpm 

Mode  1  (2.  67  deg/sec) 

J 

. 

0.  HO!) 

r  pn . 

Mod**  2  (5.  34  deg/aec* 

i  Accelerometer  set: 
j  factor 

L_ 

K 

10T 

•  onnts  /g  -  sec 
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2. 


ACCELERATION  ALGORITHMS 


Appendix  E  presents  the  development  of  the  algorithm  for  the  input 
axis  acceleration  sensed  by  the  LOGACS  accelerometer.  The  final  result  is 


A  .  =  A  cos  <t>  +  A  sin  <P  -  r 
si  xcg  t  ycg  t  a 


[<*t  +  u>z)2  +  (cos  *tu,y  -  sin  -  w2] 


X  sin  4> 

ag  [  t 


2  2  2 

(u;  +  U»  U)  )  "  COS  (w  +  U)  -  Ui  ) 

z  x  y  t  y  z  o 


cos  4>t(-u 


•  x  \  .  x  ,  2  2  2. 

w  +  u  w  )  -  sin  9  (w  r  w  -  u»  ) 


z  x  z 


t  X  z  o 


+  Z  [cos  4>  (u  u>  +  w  +  3w20)  -  sin  <f>  (ti  -  u>  w  +  3w2^) 
ag  [  t  x  z  y  o  ^t  x  y  z  o  'J 


where 

A  .  -  total  sensed  acceleration  along  accelerometer  input 
81  axis 

A  =  component  of  sensed  acceleration  along  x  due  to 
xc®  external  forces 

A  =  component  of  sensed  acceleration  along  y  due  to 
external  forces 


r  =  radius  arm  of  LOGACS 
a 

<t>t  =  turntable  rotational  rate 
<t>t  =  turntable  phase  angle 


=  inertial  rotational  rates  about  x,  y,  and  z, 

3  1/2 

=  instantaneous  orbital  rate  =  (p/R  ) 


respecti  vely 


inertial  rotational  accelerations  about  x,  y,  and  z, 
respectively 


X  ,  Y  ,  Z  =  coordinates  of  LOGACS  with  respect  to  Agena  eg 
ag  ag  ag  r  6  6 


=  (X  -  X  ),  (Y 
ac  eg  ac 


Y  ),  (Z 
eg  ac 


A  ) 
eg 


6  =  Agena  pitch  attitude  with  respect  to  local  horizontal 


=  Agena  roll  attitude  with  respect  to  local  horizontal 

Excluding  the  dynamics  of  the  LOGACS  accelerometer,  the  truncated  output 
in  units  of  counts  each  second  is 


V  (t  +  1)  =  -K 
s 


(Agi  +  A^  +  •  •  *  )  dt  +  R 


(18) 


where 

K  =  instrument  scale  factor  in  counts/g-sec 
Agj  =  sensed  acceleration  along  instrument  input  axis 
A^  =  instrument  bias 

+  •••  =  represents  nonlinear  terms,  cross-axis  acceleration 

sensitivity,  etc,,  which  are  assumed  zero  for  this  analysis 

R  =  remainder  term  of  fractional  counts  from  previous  cycle 

This  incremental  integration  is  performed  and  recorded  each  second;  there¬ 
fore,  the  number  of  counts  recorded  represents  an  average  input  g  level  over 
the  1-sec  interval.  Denoting  this  as  N,  it  is  desired  to  determine  the  accel¬ 
eration  due  to  drag  from  the  measurement(s)  N,  Using  Eq,  (17),  one  ce.n 
determine  the  drag  acceleration  over  each  1-sec  interval  from  (positive  drag 
equals  -Axcg) 
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where  all  terms  on  the  right  side  of  Eq.  (19)  are  considered  averaged  values 
over  the  1-sec  interval.  It  is  apparent  that  singularities  exist  in  this  equa¬ 
tion  at  4^  =  ±90  deg.  Since  only  N,  $t,  0,  and  ^  are  explicitly  measured,  it 
is  necessary  to  determine  the  other  terms  implicitly.  The  remainder  of  this 
section  derives  the  order  of  magnitude  of  each  term  and  presents  the  approxi¬ 
mations  made  to  these  algorithms  for  the  purpose  of  data  processing, 

3.  AGENA  RIGID  BODY  EQUATIONS  OF  MOTION 

It  is  assumed  that: 

a.  The  configuration  is  a  rigid  body  (excludes  propellant 
motions  inside  the  tanks) 

b.  The  principal  axes  are  parallel  to  the  reference  axes 

c.  The  only  external  forces  ard  moments  are  caused  by 
aerodynamics  and  control  system  gas  valve  thrusts 

d.  There  are  no  internally  generated  forces  or  moments 

e.  The  gravity  gradient  torques  are  negligible 

On  the  basis  of  these  assumptions,  the  equations  can  be  written  as  follows: 
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a.  Moments 


I  w  =  (Z  -  Z  )(T .  -  T,)  +  (Z  +  Z  )(T,  -  T,  ) 
xx  gv  eg  '  i  3  '  gv  eg  '  4  6 


+  Yeg(T2  -  T5^ 


I  W  =  (I  -  I  )w  W  +  (X  -  X  )(T,  -  Tc ) 
y  y  z  x  x  z  '  gv  eg  2  5' 


+  F  (X  -  X  )  -  F  (Z  -  Z  ) 
az  eg  cp  ax  eg  cp 


I  w  =  (I  •  I  )w  u  MX  -  X  )(T,  +  T.  -  T,  -  T, ) 
z  z  X  y  X  y  '  gv  eg  3  4  16 


F  (X  -  X  )  +  F  (Y  -  Y  ) 
ay  eg  cp  ax  eg  cp' 


where 


L  =  moments  of  inertia  about  pi  incipal  axes  (i  =  x,  y ,  z) 
w.  =  angular  accelerations  about  principal  axes  (i  =  x,  y,  z) 
w.  =  angular  rates  about  principal  axes  (i  =  x,  y,  z) 

T.  =  gas  valve  forces  (i  =  lf  2,  •••  6) 

i  =  coordinates  of  gas  valve  (i  =  X,  Yf  Z) 

8  ^ 

i  =  coordinates  of  center  of  gravity  (i  =  Xf  Y,  Z) 
cg 

i  =  coordinates  of  aerodynamic  center  (i  =  X,  Y,  Z) 
cp 

F  .  =  aerodynamic  forces  along  reference  axes  (i  =  x,  y,  z) 
ai 
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Note  that  the  moment  of  inertia  coupling  term  drops  out  of  the  equation 
of  Eqs.  (20)  because  I  =  1^.  The  aerodynamic  forces  are  determined  from 


Fax  =  -CdQS 


F 

ay 


=  cNpQsp 


F 


az 


=  -CNoQS* 


(21) 


where 

S  =  reference  area 

Cp  -■  drag  coefficient  (assumed  equal  to  axial  coefficient) 

=  normal  force  coefficients  (i  =  f>,P) 

a ,  P  =  angle  of  attack  and  side  slip  angle,  respectively 

2 

Q  =  dynamic  pressure  =  ( l / 2 ) p  V  (p  is  atmospheric  density  and 
Va  is  relative  velocity) 

b.  Sensed  eg  Accelerations 


A 


xcg 


F 

ax 

m 


F  +  (T-  +  T,  -  T, 
A  a  . aX.-  .i - 1 - L 

yeg  m 


V 


(22) 


A 

zcg 


F  +  (T,  -  T.) 
az  5  2 

m 


where  m  is  the  vehicle  mass  and  the  other  terms  are  defined  above. 
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4. 


AGENA  CONTROL  SYSTEM  RESPONSE  CHARACTERISTICS 


This  section  presents  a  tutorial  discussion  of  the  Agena  control 
system  responses  under  the  influence  of  torques,  and  develops  the  order  of 
magnitude  of  Ihe  angular  rates  and  the  angular  and  translational  accelerations 
induced  by  the  control  system.  Reference  3  presents  a  more  detailed 
description  of  the  control  system  and  its  characteristics  than  is  given  here. 

The  basic  elements  of  the  control  system  are: 

a.  Horizon  sensor  for  detecting  vehicle  pitch  and  roll 
attitude  with  respect  to  the  local  horizontal 

b.  Pitch,  yaw,  and  roll  integrating  gyros  for  detecting 
vehicle  attitude  displacements  from  the  gyros'  reference 
axes 

c.  Six  cold  gas  control  valves  for  producing  control  moments 
about  the  pitch,  yaw,  and  roll  axes  (see  Fig.  19) 

d.  Cold  gas  supply  and  associated  pressure  regulators 

e.  Associated  electronics  including  a  programmer 

The  gas  valve  with  electronics  represents  a  nonlinear  element  in  the  control 
loop  in  that  the  valve  mode  is  either  on  or  off.  Figure  20  presents  the  static 
performance  characteristic  in  which  a  given  error  voltage  (derived  from  the 
gyro  error  signals)  produces  a  given  pulsewidth  and  pulse  repetition  rate. 

The  mod  factor  is  a  product  of  the  two  curves.  It  is  seen  that  for  rather 
large  input  voltages  the  pulsewidth  remains  relatively  constant  whereas  the 
pulae  repetition  rate  increases  rapidly  with  applied  voltage.  Preceding  this 
circuit  is  a  deadband  circuit  in  which  the  voltage  must  exceed  a  given  value 
before  it  will  be  detected  by  the  control  valve  circuits.  It  should  be  noted 
that  the  minimum  valve  firing  duration  is  20  msec.  Also,  the  error  voltage 
is  a  composite  signal  derived  from  passing  the  gyro  error  signal  through  a 
lead  circuit;  therefore,  the  voltage  contains  error  rate  components  as  well. 
The  deadband  is  set  such  that  a  nominal  ±1  deg  error  at  zero  rate  will  cause 
a  minimum  firing  rate  of  the  appropriate  valve. 

The  pitch  valves  (2  and  5)  respond  to  the  pitch  gyro  error,  where  a 
pitch -up  error  will  cause  valve  2  to  fire  to  produce  a  nose -down  moment. 
Conversely,  a  pitch-down  error  will  cause  valve  5  to  fire  to  produce  a  nose- up 
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moment.  Valves  l,  3,  4,  and  6  respond  to  roll  and  yaw  gyro  error  signals 
in  the  following  manner:  for  a  yaw  right  error,  valves  3  and  4  will  attempt 
to  fire  to  produce  a  yaw  left  moment,  whereas  for  a  roll  clockwise  error, 
valves  3  and  6  will  attempt  to  fire  to  produce  a  counterclockwise  moment; 
similarly  for  a  yaw  left  error,  valves  1  and  6  will  attempt  to  fire,  whereas 
for  a  roll  counterclockwise  error,  valves  l  and  4  will  attempt  to  fire.  That 
is,  valves  1  and  3  respond  to  the  sum  of  the  roll  and  yaw  errors  whereas 
valves  4  and  6  respond  to  the  difference  of  these  errors.  As  a  result  of  this 
summing  (and  differencing),  the  probability  that  two  valves  will  fire  simul¬ 
taneously  is  very  small  because  it  occurs  only  when  one  gyro  error  (plus  rate 
component)  is  zero  when  the  other  reaches  the  voltage  level  to  cause  a 
commanded  firing.  Under  the  most  ideal  conditions  of  no  cross -axis  coupling, 
no  eg  offsets,  and  no  aerodynamic  torques,  the  pitch  moment  equations  can 
be  written  as 


I  u  =  (X  -  X  )(T  -  T.) 

y  y  gv  cgM  2  5' 


An  idealized  time  history  assuming  perfect  electronics  and  ideal  initial 
conditions  would  be  depicted  as  shown  in  Fig.  21.  During  the  interval  of  a 
gas  valve  firing  one  finds 
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For  a  control  valve  firing  duration  (At)  of  20  msec,  w  equals  0,001  deg/sec. 
Under  these  conditions  there  would  be  a  gas  valve  firing  every  2000  sec. 
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Fig.  21.  Idealized  Time  History 

When  it  is  assumed  that  there  is  a  small  constant  to'que  applied 
(maintaining  the  same  characteristic  of  minimum  pulsewidth  and  symmetrical 
rate  change  at  the  occurrence  of  the  pulse),  the  time  history  characteristics 
change  to  those  shown  in  Fig.  22. 

For  the  case  where  the  attitude  just  reaches  the  deadband,  but  does 
not  cause  a  gas  valve  firing,  the  period  between  firings  would  be  8000  sec. 
Increasing  the  applied  torque  decreases  the  period  between  firings  and 
decreases  the  attitude  excursion.  For  example,  a  torque  large  enough  to 
cause  a  firing  once  per  second  (holding  all  assumptions  including  At  =  0.02) 
would  result  in  an  attitude  excursion  of  0.  00025  deg. 

The  roll  and  yaw  responses  are  coupled  and  cannot  be  treated  on  a 
single -axis  basis,  but  the  order  of  magnitude  of  the  terms  can  be  established. 
Since  the  lever  arm  and  inertia  are  the  same  for  yaw  as  for  pitch,  the  rate 
change  and  acceleration  due  to  a  roll/yaw  valve  firing  is  the  same  for  yaw  as 
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Fig.  22.  Time  History  Assuming  a  Small  Constant  Torque 

for  pitch.  The  roll  response  for  a  roll/yaw  valve  firing  is  larger  due  to  the 
small  roll  inertia.  The  minimum  roll  acceleration  is 

.  _  0. 5  x  2.  166  ,  „  .  t/>-3  ,  /  2  „  ,,  ,  ,  2 

wx  - 172 -  =  6.  3  >  10  rad/  sec  =  0.  36  deg/sec 

and  the  corresponding  roll  rate  is 

w*  -<*>”»**  w  At  =  0.  0036  deg/ sec 

The  conditions  and  assumptions  applied  in  the  above  analysis  are  not 
generally  met  under  flight  conditions;  however,  the  results  serve  to  illustrate 
the  response  characteristics  and  to  establish  approximate  values  for  angular 
rates  and  accelerations  of  the  vehicle.  The  telemetric  data  were  used  to 
establish  the  bounds  on  the  angular  rate  terms,  and  to  verify  the  magnitude 
of  angular  accelerations  induced  by  control  gas  valve  firings. 
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Two  sources  of  telemetric  data  are  available:  that  which  was 
directly  associated  with  the  LOGACS  experiment  and  the  normal  Agena 
telemetric  data.  The  LOGACS  data,  which  were  tape  recorded  and  played 
back  during  real-time  acquisitions,  include  gyro  error  data  (pitch,  yaw,  and 
roll)  and  control  system  gas  valve  firing  data.  The  gyro  error  data  are 
recorded  at  a  2-sec  sample  interval.  For  ^ach  of  the  control  valves  the  data 
recorded  indicated  whether  one  (or  more)  firing(s)  of  the  valve  had  occurred 
in  a  1-sec  interval.  The  number  of  firings  and  the  duration  of  each  over  the 
1-sec  interval  are  not  recorded  with  the  LOGACS  data.  The  Agena  data, 
available  only  during  real-time  acquisitions,  include  a  continuous  gas  valve 
time  history  so  that  the  on-off  state  of  each  valve  can  be  deterrrured.  These 
data  can  be  used  to  establish  the  number  of  firings  and  the  time  duration  of 
each  for  the  periods  of  real-time  acquisition. 

A  rigorous  analysis  of  the  control  system  data  has  not  been  made; 
however,  from  a  cursory  examination,  various  response  characteristics 
have  been  established.  LMSC  analyzed  the  reil-time  Agena  control  valve 
data  (Ref.  4).  In  this  analysis,  data  from  Revs.  20,  25,  38,  54,  and  66 
were  used.  It  was  found  that  each  valve  firing  observed  was  of  a  minimum 
duration  (20  msec).  The  data  available  on  Rev.  66  wtre  also  analyzed 
for  pulse  repetition  rate.  Initiation  of  the  real-time  acquisition  on  this  pass 
occurred  approximately  200  sec  after  perigee.  Thus,  the  control  torques 
required  to  compensate  for  the  aerodynamic  moments  are  typical  of  the 
maximum  required  over  the  complete  mission.  The  analysis  of  this  data  set 
showed  the  average  period  between  firings  for  each  valve  (1,  2,  and  6  were 
the  only  valves  firing)  was  greater  than  1  sec.  There  was  only  one  occurrence 
in  this  data  set  where  the  period  between  firings  was  less  than  1  sec.  Based 
on  this  analysis,  it  is  assumed  for  this  report  that  all  indicated  valve  firings 
by  the  LOGACS  data  set  represent  a  single  firing  of  20  msec  duration. 

The  LOGACS  gyro  error  data  were  analyzed  in  regions  where  the 
control  lystem  was  not  operating  on  its  deadband  (see  Appendix  G,  Fig.  G-2, 
which  includes  plots  of  the  attitude  error  time  histories).  At  selected  apogee 
regions  (minimum  disturbing  torques),  the  attitude  error  rate  change  when  a 


2-14 


gas  valve  fired  was  determined.  Within  the  accuracy  of  the  telemetric  data 
and  method,  the  rate  changes  agreed  with  that  predicted  by  the  idealized 
analysis,  confirming  the  magnitudes  of  the  angular  accelerations. 

The  rate  deduced  from  analysis  of  the  attitude  error  data  reflects 
the  difference  between  the  vehicle  rate  (uu,  i  =  x,  v,  z)  and  the  gyro  refer¬ 
ence  rate.  This  rate  is  defined  as  a  limit  cycle  rate,  and  the  gyro  error 
data  were  examined  to  determine  the  maximum  limit  cycle  rate.  Relatively 
large  limit  cycle  rates  occurred  when  the  vehicle  entered  the  perigee  region. 
Here  the  build-up  of  aerodynamic  torque  drives  the  attitude  toward  the  control 
deadband.  Until  the  deadband  is  reached,  there  is  no  compensating  control 
torque.  The  maximum  rate  achieved  depends  on  the  initial  conditions  when 
entering  this  region.  Also,  for  the  later  revs.  (Revs.  47  and  up),  relatively 
large  rates  are  observed  in  the  perigee  regions  when  anomalies  occur  in  the 
yaw/roll  responses  (see  Section  D).  In  either  of  the  above  mentioned  regions, 
the  limit  cycle  rates  did  not  exceed  0.01  deg/sec  and  in  general  were  much 
smaller. 

In  the  perigee  regions  where  the  aerodynamic  torques  are  sufficient 
to  cause  the  control  system  to  operate  continuously  at  the  pitch  and  yaw  dead¬ 
bands,  vehicle  rates  cannot  be  observed  in  the  telemetric  data.  However, 
from  the  character  of  the  pulse  repetition  rate  observed  in  the  LMSC 
analysis,  it  is  concluded  that  the  system  is  operating  as  described  for  the 
torqued  ideal  system;  i.e.  ,  the  maximum  pitch  and  yaw  limit  cycle  rates  are 
0.001  deg/sec  and  the  maximum  roll  rate  is  0.  004  deg/sec. 

The  gyro  references  are  controlled  by  the  horizon  sensor  in  a  gyro- 
compassing  mode  (see  Ref,  3)  so  that  gyro  drifts  are  compensated  by  signals 
from  the  horizon  sensor  and/or  yaw  attitude  misalignment  with  respect  to  the 
orbital  plane.  The  pitch  gyro  is  programmed  at  a  constant  rate  to  compen¬ 
sate  for  the  orbital  rate.  Differences  between  the  instantaneous  orbital  rate 
and  the  programmed  rate  are  also  compensated  by  the  horizon  sensor.  Thus, 
the  limit  cycle  rates  deduced  from  telemetry  for  roll  and  yaw  are,  to  first 
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order,  the  vehicle  rates  w  ,  u>  .  Also,  to  first  order,  the  pitch  vehicle  rate 

X  z 

is  the  sum  of  the  orbital  rate  and  the  limit  cycle  rate,  i.e. 


0 J 


w  +  6 
o 


The  translational  acceleration  during  a  gas  valve  firing  is 


A.  =  — j- =  =  0.  0835  ft/sec2 

leg  m  59. 5 

Summarizing,  the  roll  and  yaw  vehicle  rates  do  not  exceed  0.  01  deg/ 
sec  and  the  vehicle  pitch  rate  is  bounded  within  the  limits  of  orbital  rate  plus 
or  minus  0.01  deg/ sec.  The  accelerations  experienced  during  a  gas  valve 
firing  are  summarized  in  Table  10. 

Table  10.  Accelerations  During  Gas  Valve  Firing 
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5. 


AERODYNAMIC  FORCES  AND  MOMENTS 


Reference  5  presents  aerodynamic  data  for  the  Agena  configuration 
flown  on  this  mission.  In  preparing  the  data,  a  flat-plate  front  end  con¬ 
figuration  was  assumed.  The  data  are  presented  in  terms  of  lift,  drag,  and 
moment  coefficients.  For  this  analysis  it  is  more  convenient  to  deal  with 
body  axis  coefficients  (normal  and  axial  force)  and  a  center  of  pressure 
distance  fot  determining  moments.  The  equations  used  for  this  transforma¬ 
tion  are: 


Cj^  =  C^  cos  a  +  Cp  sin  a 


C  .  =  C„  cos  o  -  G.  sin  o  =  C„ 
AD  L  D 


L  (-C  /C„.)L 

cp  '  m  N  r 


(normal  force  coefficient) 

(axial  force  coefficient) 

(center-of-pressure  distance 
from  Station  247) 


=  c  N<«)/» 


(normal  force  coefficient  slope) 


The  data  of  Ref.  5  for  an  exospheric  temperature  of  1500°K  were  used  to 
prepare  the  set  of  aerodynamic  terms  shown  in  Table  11. 


Table  11.  Aerodynamic  Coefficients 


h,  nmi 


CD  {<>  =  0) 

CNo  (per  deg) 

L  ,  ft 
cp’ 

X  ,  i 
CP 

2.  3 

0. 0275 

-  6.7 

2.0 

2.56 

0.  048 

-  7.6 

1.  1 

2.68 

G,  052 

-  9.4 

-0.7 

2.83 

0.  056 

-11.5 

-2.8 

2.  95 

0.  059 

-U.5 

-2.8 

Above  95  nmi,  the  x  axis  location  of  the  aerodynamic  center  of  pressure  is 
predicted  to  be  6ehind  the  eg,  providing  aerodynamic  stability.  The  flat- 
plate  front  end  configuration  places  the  aerodynamic  center  for  axial  force 
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on  the  x  axis.  Lateral  eg  offset  provides  a  lever  arm  for  the  axial  force, 
resulting  in  a  moment.  The  combination  of  off-axis  eg  and  stable  aero¬ 
dynamic  normal  forces  results  in  a  shift  of  the  trim  angle  of  attack  from 
zero  to 


a 


C.Z 

_A_££. 

CL.  X 
No-  cp 


provided  X  is  negative.  For  the  aerodynamic  and  eg  data  given,  the  trim 
cp 

angle  of  attack  exceeds  the  control  system  deadband  for  all  altitudes  so  that 
a  control  torque  is  necessary  to  compensate  the  total  induced  aerodynamic 
moment. 

Using  the  aerodynamic  coefficients  of  Table  11  and  density  data  from 
the  1962  standard  atmosphere  tables,  we  calculated  the  aerodynamic  forces 
and  accelerations  predicted  as  a  function  of  altitude;  they  are  presented  in 
Table  12.  The  density  data  and  the  drag  acceleration  in  units  of  accelerom¬ 
eter  counts  are  plotted  in  Fig.  23. 

6.  ORDER  OF  MAGNITUDE  ESTIMATES 


Based  on  the  data  presented  in  Sections  4  and  5,  an  order  of  magni¬ 
tude  estimate  can  be  rswle  for  the  individual  terms  of  the  sensed  acceleration 
given  in  Eq.  (17i.  Before  a  term-by-term  estimate  is  made,  certain  terms 
are  eliminated  from  Eq.  (17)  based  on  their  order  of  magnitude  as  follows: 


Limit  cycle  rate  terms:  Using  the  maximum  value  of 
0.01  deg/sec  (or  2  x  10*4  rad/sec),  products  of  limit  cycle 
rates  L-iuj,  u-j©;  i,  j  =  x,  z)  are  less  than  4  X  10*®  rad/sec2. 

With  lever  arms  less  than  10  ft,  the  resulting  accelerations 
arc  less  than  4  X  10'?  ft/sec2  (1.25  x  10-8  g  =  0.  125  accelerom¬ 
eter  counts).  Thus,  products  of  limit  cycle  rates  are 
eliminated  from  further  consideration. 

Orbital  rate  terms:  Using  a  maximum  value  for  the  orbital 
rate  term  (wQ)  of  4.  35  deg/min  (1.  3  X  10-3  rad/sec),  products 
of  orbital  rate  and  limit  cycle  rates  are  less  than  2,  5  X  10*^ 
rad/sec2.  With  lever  arms  less  than  10  ft,  the  resulting 
accelerations  are  less  than  2.  5  X  10"®  ft/sec^  (8  X  10'®  g  =  0,  8 
accelerometer  counts).  The  square  of  the  orbital  rate  term 
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Table  12.  Aerodynamic  Forces  and  Accelerations 


ACCELERATION,  counts 


timers  the  turntable  radius  |ra),  or  the  y-axis 

lever  arm  (Ya„),  is  less  than  6  X  10~7  ft/sec-  ?  ? 

(2  X  10'^  g  =  6.2  accelerometer  counts).  The  term  (u.y  -  u^) 
reduces  to  2->c  6  4  0^),  and  by  (a)  above  reduces  further  to 

2wcfi.  This  term  has  a  maximum  magnitude  of  ±1.6  accel¬ 
erometer  counts.  None  of  the  above  terms  is  considered 
large  enough  for  further  consideration. 

Coupling  of  yaw  rate  and  turntable  rate:  The  coupling  term 
(2ra<t>t^z'  can  result  in  an  acceleration  of  9.2  X10~6  ft/sec^ 

(3  X  10-”  g  =  3  accelerometer  counts).  Although  this  magni¬ 
tude  is  detectable,  it  occurs  only  for  mode  2  of  the  calibration 
cycle  when  the  yaw  rate  is  maximum.  Both  the  probability  of 
occurrence  and  the  magnitude  are  small  so  that  it  too  is 
eliminated  from  further  consideration. 


Gravity  gradient  terms:  The  terms  associated  with  attitude 
(3Zagw^0  and  3Zag'-5tt>)  are  of  order  ±1.69  *  10-7  ft/sec^ 
(0.05  accelerometer  counts)  and  are  eliminated  from  further 
consideration. 


As  a  result  of  the  above  considerations,  all  vehicle  ”ate  terms  in  Eq.  (17) 
can  be  eliminated  and  the  sensed  acceleration  equation  reduces  to 


A  .  =  -r  4>?  +  cos  6, (A  -  Y  u  +  Z  w  ) 
si  at  txcg  ag  z  ag  y 


t  sin  $  (A  4  X  w  -  Z  w  ) 
t  yeg  ag  z  ag  x 


Based  on  an  accelerometer  scale  factor  of  10  counts/g- sec,  and 
with  contributions  of  aerodynamic  moments  neglected  (see  Section  7),  the 
order  of  magnitude  estimates  for  the  terms  in  Eq.  (23)  are  presented  in 
Table  13. 

7.  TIME  AVERAGE  RESPONSES 

As  discussed  in  i  ection  4,  the  control  system  responds  to  limit  the 
amplitude  of  attitude  excursion.  Under  the  relatively  high  aerodynamic  torque 
experienced  below  130  nrni  altitude,  the  attitude  error  is  nearly  constant 
(well  within  0,  05  deg  of  its  deadband  limit),  and  the  limit  cycle  rates  (integral 
of  -  )  are  bounded  at  approximately  their  threshold  values.  Thus,  the  time 
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Table  13.  Order  of  Magnitude  Estimates 


WTOrf;rr,r  ri**~fW*'?7WW 


average  of  3  has  a  mean  (expected)  value  of  zero.  The  integral  of  r  x  w 
terms  (components  of  which  are  sensed  by  the  LOGACS  accelerometer)  would 
also  have  zero  mean.  In  essence,  the  mean  of  the  gas  control  torque  is  equal 
in  magnitude  and  opposite  in  polarity  to  the  mean  of  the  aerodynamic  torque. 

If  the  pneumatic  control  system  was  a  linear  system  (proportional 
instead  of  pulsed  valves),  the  r  x  u  terms  would  be  negligible.  With  the 
pulse  system,  however,  the  accelerometer  output  at  each  second  will  be 
related  to  the  difference  of  the  torque  sources  over  the  interval.  Over  a 
given  interval  the  probability  is  negligible  tl  at  the  integrated  magnitudes  of 
each  are  exactly  equal;  therefore,  there  is  a  variance  associated  with  each 
accelerometer  output.  This  variance  is  also  associated  with  the  magnitude 
of  the  aerodynamic  torque  as  well  as  the  orientation  phase  angle  during  the 
rotating  calibration  modes.  The  accelerometer's  dynamic  response  to  pulsed 
accelerations  also  contributes  to  the  variance  of  the  accelerometer's  output. 
Consequently,  a  given  numerical  value  cannot  be  established  for  the  vari¬ 
ance;  however,  Table  13  does  give  the  order  of  magnitude  of  a  maximum 
error  value.  Filtering  techniques  are  employed  to  minimize  the  variance 
associated  with  the  processed  output  data. 

B.  LOGACS  DATA  PROCESSING  METHODS 

The  telemetered  functions  which  were  processed  and  their  results 
included  in  this  report  are: 


a.  Accelerometer  counts  (N) 

b.  Control  valve  firing  occurrences  (T^) 

c.  Pitch,  yaw,  and  roll  gyro  error  (0^,  ,  4>e) 

d.  Pitch  and  roll  horizon  sensor  (6,  ) 


The  data  rate  for  the  first  two  was  once  per  second  while  the  latter  two 
had  a  data  rate  of  once  every  other  second.  The  time  base  used  for 
processing  and  presenting  the  data  is  the  LOGACS  clock  data  designated 
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TRG.  For  purposes  of  this  report  the  relationship  between  system  time 
(SYST)  and  TRG  is 

SYST  =  67,  507  -  m(86,  400)  +  TRG 

where 

m  =  0,  1,  2,  3,  4  (m  =  0  on  22  May) 

A  complete  list  of  the  telemetered  functions  and  also  the  accuracy  of  the 
TRG  clock  is  given  in  Section  I. 

There  were  a  number  of  obvious  errors  in  the  accelerometer  data 
displayed  as  a  512  or  4096  count  deviation  from  the  mean  of  the  surrounding 
data  points  due  to  transmission,  receiving,  or  decoding  problems.  Most 
probably  there  are  an  equal  number  of  1,  8,  and  64  (and  multiples  of)  count 
errors  which  are  not  recognizable.  The  larger  errors  were  edited  from  the 
source  data  tape.  The  frequency  of  occurrence  of  this  type  error  is  such 
that  the  smaller  errors  would  not  contribute  significantly  Co  the  final  results. 

Figure  F-l  of  Appendix  F  contains  an  example  of  the  raw  data  plots 
(excluding  horizon  sensor  data)  covering  the  data  span  from  280,  000  to 
286,  300  TRG  sec.  This  time  span  includes  the  apogee  regions  of  Revs. 

52  and  5  3  and  the  perigee  of  Rev.  5  3.  The  control  valve  data  are  dis¬ 
played  by  an  X  in  the  appropriate  row  fo»  each  valve  when  the  data  indicated 
the  occurrence  of  a  valve  firing.  The  accelerometer  and  gyro  functions  are 
plotted  at  their  respective  amplitudes  and  data  rates.  The  initiation  of  each 
accelerometer  calibration  cycle  mode  is  indicated  on  the  plot.  Gaps  in  the 
data  reflect  missing  data  points  either  in  data  transmission,  source  proces¬ 
sing,  or  from  the  editing  procedure  described  above. 

There  were  a  number  of  approaches  investigated  for  processing  the 
accelerometer  data  to  meet  the  objectives  of  accelerometer  calibration  and 
determination  of  drag  acceleration.  The  methods  adopted  for  this  report 
were  in  the  interest  of  expediency  with  some  loss  in  accuracy.  Three  sepa¬ 
rate  methods  were  used  which  are  described  below. 
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1.  ACCELEROMETER  CALIBRATIONS,  APOGEE 

The  calibration  of  the  accelerometer's  scale  factor  and  bias  was 
performed  using  modes  1  and  2  data  from  apogee  regions  only.  This  was 
done  because  a  minimum  of  external  disturbing  accelerations  occur  in  this 
region.  The  data  set  used  for  these  calibrations  were  edited  for  observed 
responses  due  to  gas  valve  firings  and  any  accelerometer  count  anomalies 
that  exceeded  50  counts  from  the  surrounding  data  mean.  The  editing  of  data 
due  to  gas  valve  firings  allowed  the  elimination  of  from  Eq.  (23)  as  well 
as  the  contribution  of  gas  valve  firings  to  A  .  Also,  the  aerodynamic 
normal  forces  are  negligible  in  the  apogee  region  so  that  under  these  condi¬ 
tions  the  sensed  acceleration  equation  reduces  to 


A  .  =  -r  A  +  A  cos  <i> 
si  at  xcg  xt 


The  corresponding  accelerometer  output  each  second  is 


N.  =  K 
i 


+  a  ,  cos  4>  }  dt  +  B 


.)■ 


(24) 


(25) 


With  the  initiation  of  the  first  calibration  cycle  known,  4>t  and  4>t  are 
deterministic  functions  of  time.  The  value  for  the  turntable  radius  arm  r 

a 

is  that  determined  from  ground  calibrations.  Therefore,  the  unknowns  in  the 
above  equation  are  the  scale  factor  K,  the  accelerometer  bias  B,  and  the  drag 
acceleration  A^..  To  determine  these  quantities,  we  assumed  that  the  drag 
acceleration  behaved  as  a  quadratic  in  time  over  the  mode  1  and  2  calibration 
cycle  interval  (1024  sec).  Using  this  approximation  for  drag  results  in  a  set 
of  overdetermined  linear  equations  for  the  unknown  parameters.  These  equa¬ 
tions  were  mechanized  in  a  computer  program  and  the  solutions  for  the 
unknowns  were  obtained  using  a  least  squares  formulation. 
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Having  determined  the  scale  factor  and  bias,  one  can  compute  the 
measured  drag  acceleration  from 


A,  =  — 
d  cos 


1 -  (-1 

s  4>t  \ 


N  -  B 


•  ra4>t 


Since  this  equation  has  a  singularity  at  4>^  =  ±90  deg,  calculations  using 
Eq.  (26)  were  not  performed  for  phase  angles  =  ±  ; 9 0  ±  12  degj.  Fig¬ 
ure  F-2  of  Appendix  F  presents  the  results  from  the  Rev.  52  apogee 
calibration.  Both  the  quadratic -derived  acceleration  profile  and  the  one 
calculated  using  Eq.  (26)  are  plotted.  In  the  interval  of  ±12  deg  from  the 
singularity,  A(j  is  plotted  as  a  zero  value.  This  plot  demonstrates  the 
residual  noise  on  the  raw  data.  It  should  be  noted  that  in  this  presentation 
form,  the  noise  is  amplified  by  the  factor  1/cos  $  . 

A  total  of  59  data  sets  were  processed  using  the  procedure  described 
above.  The  scale  factor  and  bias  values  determined  from  these  calibrations 
are  presented  in  Section  I.  The  acceleration  profiles  derived  from  the  poly¬ 
nomial  are  presented  in  Section  C. 

2.  ACCELEROMETER  ROTATING  MODES,  PERIGEE 

In  the  perigee  regions,  at  least  one  control  valve  was  firing  each 
second;  therefore,  the  transients  induced  by  the  valve  firings  cannot  be  edited 
without  deleting  the  complete  data  set.  As  discussed  in  Section  7,  the  angular 
accelerations  induced  by  aerodynamics  and  control  valve  firings  average  to 
zero  in  the  perigee  region.  Deleting  the  angular  acceleration  terms  (w.) 
from  Eq.  (23)  results  in  the  following  sensed  acceleration  equation 


A  .  =  -r  (  A  cos  <£.  4  A  sin  4> 
si  at  xcg  t  ycg  t 


The  corresponding  accelerometer  output  each  second  Is 


,-ttAt  /  2 

N.  =  K  /  I  r  ir"  +  A ,  cos  <>  -  A  s 

i  Jt  \  a  t  d  ycg 


u.  *t) 


dt  +  B 
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Since  estimates  of  K  and  B  are  known  from  the  apogee  calibration,  the 
unknowns  in  Eq.  (28)  are  the  drag  and  normal  acceleration  components. 
Various  approaches  can  be  used  to  determine  these  by  appropriate  modeling 
of  the  acceleration  components  and  least  squares  fitting  of  the  data. 

Since  the  primary  objective  was  the  determination  of  the  drag  accel¬ 
eration  component,  it  was  decided  to  solve  Eq.  (28)  for  the  drag  component 
assuming  the  normal  acceleration  component  was  zero.  With  this  assump¬ 
tion,  the  drag  acceleration  component  becomes  Eq.  (26)  with  an  error  term 
of  A^  tan  4>t<  The  error  term  is  zero  or  negligibly  small  in  the  immediate 
region  of  =  0  or  180  deg.  All  the  raw  data  were  processed  using  Eq.  (26) 
and  plotted.  For  the  mode  1  and  2  data,  an  inverted  triangular  symbol  was 
placed  on  the  plot  when  the  phase  angle  was  0  or  180  deg.  A  typical  plot  is 
shown  in  Fig.  F-3  of  Appendix  F  which  covers  the  same  time  interval  as  die 
raw  data  set  of  Fig.  F-l.  The  data  are  corrected  for  bias,  and  the  amplitude 
is  normalized  in  units  of  accelerometer  counts  using  a  normalizing  scale 

7 

factor  cf  10  counts /g-sec.  These  plots  were  then  used  to  hand-reduce  the 
perigee  rotating  mode  data.  The  eyeball  average  around  the  indicated  phase 
angles  of  0  or  180  deg  was  recorded  and  the  results  of  this  processing  are 
presented  in  Section  C.  These  plots  were  also  useful  for  processing  the 
fixed  mode  (3  and  4)  data,  which  is  described  in  Section  3. 

For  the  data  set  presented  in  Fig.  F-3,  the  mode  1  and  2  results 
appear  well  behaved  (282,  240  <  TRG  <  283,  270  sec).  This  is  due  to  a 
response  anomaly  experienced  during  this  time  interval  in  which  the  normal 
acceleration  component  (A^c^)  was  near  zero  (see  Section  D).  Figure  F-4  of 
Appendix  F  presents  a  more  typical  time  history  plot  of  the  perigee  mode  1 
and  2  corrected  data.  The  influence  of  the  normal  acceleration  component  is 
immediately  apparent  as  well  as  the  difficulty  in  extracting  the  drag  accel¬ 
eration  component. 

3.  FIXED  MODE  DATA  PROCESSING 

For  the  fixed  orientation  modes,  Eq.  (26)  also  holds  where 
* 

cos  -  ±1  and  ^  =  0.  A  linear  least  squares  filter  was  used  tc  process 
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these  data  sets.  The  filter  time  interval  was  chosen  as  20  sec  in  the  regions 
where  the  drag  acceleration  exceeded  60  accelerometer  counts,  and  60  sec 
when  it  was  less  than  60  counts.  The  data  plots  of  corrected  acceleration 
were  used  to  set  up  the  processing  intervals  and  to  edit  data  points  or  inter¬ 
vals  which  exhibited  questionable  data. 

One  notable  data  time  interval  when  the  accelerometer  readings  did 
not  indicate  the  acceleration  level  was  the  interval  immediately  following  the 
mode  transitions  from  rotating  to  fixed  orientation  (Modes  2  to  3  and  3  to  4) 
in  the  apogee  regions.  The  time  required  to  traverse  the  instrument  dead 
zone  for  a  change  in  polarity  of  the  total  force  on  the  proof  mass  (accelera¬ 
tion  plus  bias)  is  inversely  proportional  to  the  square  root  of  the  force  level. 
Knowing  the  bias  term,  we  could  calculate  the  acceleration  from  the  time 
interval.  For  this  report,  these  intervals  were  edited  from  the  data 
processing. 

The  results  from  the  fixed  mode  processing  are  also  presented  in 
Section  C. 

4.  CONTROL  SYSTEM  FUNCTIONS 

The  control  functions  of  gyro  error  data  (pitch,  yaw,  and  roll)  and 
horizon  sensor  data  (pitch  and  roll)  were  processed  using  a  linear  least 
squares  filter  and  a  30-sec  processing  interval.  No  attempt  was  made  to 
edit  extraneous  data  points.  The  differences  between  the  filtered  pitch  gyro 
and  pitch  horizon  sensor  data,  and  the  corresponding  roll  differences,  were 
also  determined  and  plotted. 

The  control  valve  firing  data  were  processed  making  the  assumption 
that  each  indicated  valve  firing  was  representative  of  a  single  firing  for  a 
minimum  20-msec  duration  (see  Section  4).  The  accumulated  number  of 
firings  over  a  100-sec  time  interval  was  determined  for  each  valve  (N.). 
Also,  the  data  duty  cycle,  i.  e.,  the  number  of  recorded  data  times  (M)  in 
the  100-sec  interval,  was  determined.  Using  the  predicted  values  for  the 
gas  valve  thruet(s)  level  and  the  vehicle  mass,  we  determined  the  average 


accelerations  due  to  control  firings  along  the  vehicle  y  and  a  axes  in  units 
of  accelerometer  counts  as  follows 


A  =  K 
ycg  c 


[N3tfV  (Nj  run 


N5-N2 


where 


Kc  =  K(^)  =  ‘»7(H^)  =  «-3 


The  time  associated  with  each  processed  interval  was  computed  as 
the  time  average  for  the  data  set,  i.e., 


\'TRG. 

TRG  =  —  , 
a  M 


Since  the  processed  valve  firing  counts  N.  cannot  exceed  M,  the  above 
relationships  have  a  saturation  level  of  *52.  3  for  A  and  ±104.  6  for  A 

The  average  roll  control  moment  was  calculated  based  on  the  pre¬ 
dicted  vehicle  eg  data  and  the  assumed  value  for  control  valve  thrust(s)  as 
follows 

Mr  *  W  [,Zgv  •  ^cg1  ,N1  -  N3>  *  Zgv  +  zcg»  ,N4  -  N6> 

+  ,Ycg>  <NZ  -  N5>] 


_  0.01 


i \a  [2*041  -  N3>  t  2.291  (N4  -  N6)  +  0.1915  (N 


^2  *  n5)  (30) 
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A  bias  in  the  reduced  data  would  indicate  an  error  in  the  assumed 
eg  data.  Systematic  trends  in  the  perigee  region  would  be  indicative  of  aero¬ 
dynamic  roll  moments. 

The  results  of  the  control  system  data  processing  are  presented  in 
Section  C. 

C.  DATA  PROCESSING  RESULTS 

The  drag  acceleration  profile  reduced  from  the  data  is  plotted  in 
Fig.  G-l  of  Appendix  G.  The  acceleration  is  normalized  into  units  of  counts 
where  one  count  represents  10  ^  g.  A  different  symbol  is  used  dependent  on 
the  processing  method  (see  Section  B)  as  follows: 

Method  Plot  Symbol 

Perigee  rotating  Modes  1  fit  2  A 

Fixed  Modes  3  &  4  + 

Apogee  rotating  Modes  1  &  2  V 

Included  in  this  plot  are  data  from  the  best -fit  ephemeris  (Ref.  2)  of  latitude, 
longitude,  and  altitude,  and  the  drag  acceleration  determined  from  the  model 
atmosphere  used. 

There  are  some  redundant  regions  of  processed  data  where  the  hand- 
reduction  method  was  used  to  read  data  from  the  fixed  orientation  modes  and 
from  the  apogee  calibration  modes.  The  gaps  in  the  data  reflect  data  voids. 
The  discontinuities  of  ephemeris  data  reflect  the  differences  between  the 
independent  vector  solutions  (see  Ref.  2). 

The  data  associated  with  the  control  system  are  plotted  in  Fig.  G-2 
of  Appendix  G.  The  drag  acceleration  is  also  plotted  so  that  the  charac¬ 
teristics  of  response  anomalies  can  be  observed  more  easily.  Plotted  are 
the  data  duty  cycle,  gyro  errors  (pitch,  yaw,  and  roll),  differences  between 
gyro  error  and  horizon  sensor  (pitch  and  roll),  roll  moment,  negative  z  and  y 
axis  acceleration  due  to  control  valve  firings,  and  the  drag  acceleration 
profile. 
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Wit!1,  respect  to  the  results  obtained  for  the  normal  acceleration 
components,  nowhere  in  the  processed  data  set  did  the  results  produce  a 
saturated  value  for  these  components.  This  result  lends  credence  to  the 
premise  stated  in  Section  4  that  the  minimum  average  control  valve  firing 
period  was  greater  than  1  sec. 

The  accuracy  of  the  derived  normal  accelerations  degrades  with 
reduced  data  duty  cycle.  For  example,  if  the  duty  cycle  were  1  percent 
[M  =  1  in  Eq,  (29)]  for  a  100- sec  processing  interval,  then  the  indicated 
occurrence  (or  nonoccurrence)  of  a  firing  would  produce  a  saturated  (or  null) 
result.  For  this  example,  the  probability  that  the  result  is  correct  (assum 
ing  no  other  contributing  factors)  is  0.01,  since  nothing  is  known  about  the 
control  valve  firing  state  for  the  other  99  sec.  For  duty  cycles  between  1 
and  100  percent,  the  probability  that  the  reduced  data  represent  an  accurate 
estimate  is  also  dependent  on  the  duty  cycle  data  sequence.  The  associated 
probabilities  with  a  given  duty  cycle  have  not  been  calculated  here,  but  this 
accuracy  degradation  is  pointed  out  so  that  caution  is  exercised  about  inter¬ 
pretations  placed  on  the  reduced  normal  acceleration  data  for  regions  of  low 
data  duty  cycle.  To  a  lesser  degree,  low  duty  cycles  also  affect  the  results 
for  the  accelerometer  data,  particularly  if  the  cause  of  the  low  duty  cycle 
was  a  low  telemetering  signal  strength. 

The  accumulated  number  of  valve  firings  for  the  data  set  processed 
is 

Nj  =  18,684 

N2  =  21,648 

N3  =  620 

N,  -  145 
4 

N5  =  79 

N6  =  16,469 

Nt  =  57,645 
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The  accumulated  data  duty  cycle  is  63.666  percent.  Assuming  the 
characteristics  of  valve  firings  were  the  same  for  the  time  intervals  not 
recorded,  the  estimated  total  number  of  firings  for  the  mission  from  injec¬ 
tion  is  90,  500.  Assuming  a  specific  impulse  of  50  sec  for  the  gas  mixture 
used,  the  control  gas  weight  expended  was 

TAt 

w  =  nt  r1  =  18  lb 

sp 

D.  RESPONSE  ANOMALIES 

The  data  processing  performed  for  this  report  was  intended  to 
establish  gross  trends  from  which  regions  of  particular  interest  could  be 
isolated  for  detailed  analysis.  One  region  that  exhibits  a  pronounced 
anomaly,  observed  as  a  large  reduction  in  drag  from  that  expected  on 
Rev.  47  and  following,  is  when  the  vehicle  is  approaching  perigee.  This 
phenomenon  is  most  clearly  displayed  in  the  perigee  data  of  Rev.  53  [see 
Fig.  G-l(c)].  The  raw  data  for  this  pass  are  included  in  Fig.  F-l  and  the 
corrected  accelerometer  data  are  presented  in  Fig.  F-3. 

Correlated  with  the  reduction  in  drag  acceleration  is  a  change  in 
polarity  of  yaw  control  force  (see  Fig.  G-2),  and  consequently  a  correspond¬ 
ing  polarity  change  of  yaw  control  moment.  From  this  response  charac¬ 
teristic,  and  the  assumption  that  it  is  externally  caused,  it  must  be  concluded 
that  during  this  interval  the  yaw  aerodynamic  moment  had  changed  polarity. 
The  equation  for  yaw  aerodynamic  moment  is  [see  Section  3,  Eq.  (20)] 
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Since  the  second  term  in  the  moment  equation  did  not  change  sign,  the 

change  in  polarity  must  be  attributed  to  the  aerodynamic  side  force  term  (F  ). 

ay 

This  in  turn  leads  to  the  conclusion  that  there  occurred  a  large  change  in  the 
yaw  angle  of  attack  ( p )  caused  by  a  large  horizontal  wind  component.  A 
similar  phenomenon  of  lesser  degree  is  observed  in  the  pitch  axis  immedi¬ 
ately  following  perigee  of  Rev.  53.  By  the  same  analogies,  it  is  con¬ 
jectured  that  a  vertical  component  of  wind  was  encountered.  Accepting  the 
existance  of  large  wind  components,  the  effect  of  winds  on  the  dynamic  pres¬ 
sure  term  (Q)  must  also  be  investigated  to  fully  evaluate  the  drag  acceleration 
reduction  and  the  inferred  atmospheric  density. 

Before  the  anomalous  regions  can  be  analyzed  comprehensively,  a  more 
precise  data  reduction  procedure  will  be  required  in  which  the  total  lateral 
acceleration  is  determined  from  the  calibration  rotating  mode  data.  Accom¬ 
panying  this  is  the  requirement  for  a  more  precise  determination  of  aero¬ 
dynamic  coefficients  in  which  the  open  front  end  configuration  is  used  as  the 
model  instead  of  the  assumption  of  a  flat  plate.  Also,  calculations  for  larger 
angles  of  attack  than  considered  for  the  data  of  Ref.  5  will  be  required.  With 
this  information  and  assuming  aerodynamic  symmetry  in  pitch  and  yaw,  the 
magnitude  of  Q  and  the  vertical  and  horizontal  components  of  V&  can  be 
determined  (at  least  in  principle).  The  remaining  unknowns  then  relate  to 
density  (p)  and  the  in-track  component  of  V  .  Unfortunately,  there  was  no 
vehicleborne  instrumentation  which  could  be  used  to  determine  either:  how¬ 
ever,  a  bounded  functional  relationship  is  established.  With  reduced  time 
histories  of  the  functional  relationship  in  conjunction  with  analysis  of 
atmospheric  dynamics,  the  atmospheric  properties  sought  can  be 
determined. 

Another  response  anomaly  observed  in  the  data  is  associated  with 
the  reduced  difference  between  the  roll  gyro  and  roll  horizon  sen sor  data 
(see  Fig.  G-2).  A  cursory  examination  leads  to  the  conclusion  that  ,  scale 
factor  of  approximately  20  percent  exists,  either  in  the  control  system  or  the 
telemetry  of  the  respective  functions.  To  isolate  the  source  would  require  a 
dynamic  analysis  of  the  gyrocompassing  loop  under  the  influence  of  the 
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observed  response.  The  utility  of  this  function  is  to  determine  yaw  gyro 
drift  and  its  effect  on  the  yaw  reference  of  the  control  system  (see  Ref.  3 
for  a  description  and  analysis  of  the  gyrocompass  loop  and  sources  of  yaw 
reference  error).  The  anomaly  is  not  of  sufficient  magnitude  to  cause  great 
concern,  and  the  function  magnitude  implies  the  yaw  gyro  drift  was  well 
within  its  tolerance  bounds. 
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APPENDIX  A 


EXPERIMENT  CLOCK  ERROR 


The  experiment  clock  was  a  model  11477  Time  Reference  Generator 
made  by  Adcole  Corporation.  The  clock  played  a  key  role  in  the  LOGACS 
experiment  because  the  centrifuge  mode  change  pulses  and  the  centrifuge 
motor  supply  frequency  were  derived  from  the  clock  circuitry. 

After  the  LOGACS  flight,  the  local  time  of  eac  station  acquisition 
was  correlated  with  the  experiment  clock  time  acquired  during  each  station 
acquisition.  A  plot  of  the  relationship  between  the  local  station  time  and 
experiment  clock  time  is  shown  in  Fig.  A.  1.  As  can  be  quickly  seen  from 
the  figure,  the  experiment  clock  steadily  gained  time  (i.e.  ,  was  fast)  by 
about  one  part  in  50,000. 

The  clock  accuracy  as  a  function  of  temperature  was  measured 
during  the  acceptance  testing.  The  experiment  clock  was  about  65°  F  at  the 
time  of  turn-on  of  the  clock  and  then  slowly  approached  a  temperature  of 
50  ±  5°F.  The  experiment  clock  error  as  a  function  of  temperature  is  shown 
in  Fig.  A.  2.  As  can  be  seen,  the  clock  error  over  the  45°F  to  75°r  region 
is  almost  exactly  one  part  in  50,000  fast  (i.  e.  ,  the  experiment  clock  was 
indicating  a  larger  value  of  time  than  had  actually  transpired). 

The  correlation  between  the  laboratory  tests  on  the  experiment 
clock  and  the  data  obtained  from  the  clock  while  on-orbit  certainly  provide 
graphic  evidence  of  the  stability  of  the  time  reference  during  the  100  hr  of 
orbital  operation. 

Since  all  data  processing  was  accomplished  with  respect  to  the 
experiment  clock  time,  the  clock  error  did  not  really  propagate  as  an  error 
in  the  calculations  that  were  made.  The  known  clock  error  results  ivi  a  scale 
factor  error  of  one  part  in  25,000,  or  0.004  percent,  which  is  significantly 
below  the  LOGACS  experimental  error.  It  is  only  when  one  desires  to 
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correlate  the  drag  information  derived  from  the  LOGACS  data  with  the 

ephemens  of  the  satellite  vehicle  that  one  must  reckon  with  the  clock  error 

af  alU  F°r  most  such  Purposes,  the  known  clock  error  of  one  part  in  50,000 
is  small  enough  to  ignore. 
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Fig.  A.  2.  Experiment  Clock  Error  as  a  Function  of  Temperature 
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APPENDIX  C 


SPECIFICATION  FOR  LOGACS  DATA  ANALYSIS  PROGRAM 

C.  1  INTRODUCTION 

The  purpose  of  the  LOGACS  experiment  is  to  calibrate  a  low-g 
accelerometer  aboard  a  satellite  vehicle  and  use  it  to  measure  the  accelera¬ 
tion  of  the  vehicle  due  to  nongravitational  forces,  primarily  aerodynamic 
drag. 

The  accelerometer  is  calibrated  by  rotating  it  at  two  different  known 
angular  rates,  thereby  inducing  known  centripetal  accelerations.  The 
accelerometer  is  also  stopped  at  two  fixed  points  so  that  its  input  axis  is 
parallel  to  the  satellite  longitudinal  axis  pointing  either  forward  or  aft. 

The  input  to  this  data  analysis  program  is  a  floating  binary  digital 
tape.  The  output  will  be  values  of  accelerometer  scale  factor  and  bias  and  a 
time  history  of  the  vehicle  acceleration  due  to  drag. 

C.  2  SYMBOLS 

i  experiment  clock  time  expressed  as  an  integral 

number  of  seconds 

tg  mode  starting  time  (sec) 

tr  ( j )  corrected  total  of  accelerometer  pulses  accumulated 

over  20R  sec  from  tg  +  300  +  20(j-10)  to  to  +  300 
t  20(j-9)R  or  over  20  sec  from  to  +  20 (j  +  41  to  to  + 

20  ( j  +  5 ) 

ir'(i)  corrected  total  of  accelerometer  pulses  accumulated 

over  4  sec  from  i  -  2  to  i  +  2 

flr”{i)  corrected  total  of  accelerometer  pulses  accumulated 

over  1  sec  from  i  -  1  to  i 

f(i)  original  total  of  accelerometer  pulses  accumulated 

over  1  sec  from  i  -  1  to  i 

gfj)  weighting  factor  applied  toor(j)  for  interpolation 

R  ratio  between  nominal  and  actual  turntable  rates 

(R  =  184,798, 313/ 164, 450, 000) 


m(i)  calibration  cycle  mode  number  (m  =  1,2,3,  4,  or  5) 

S(m,  l)  weighted  summation  of  <r(j)  for  mode  m 

A  centripetal  acceleration  induced  by  turntable  in 

mode  l  (ft/sec2) 

f  accelerometer  scale  factor  (pulses/ft/sec) 

b  accelerometer  bias  (ft/sec^) 

d(j)  computed  values  of  drag  acceleration  (ft/sec^) 

2 

D(m)  average  value  of  drag  acceleration  (ft/sec  ) 

.  3 

D(m)  average  rate  of  change  of  drag  acceleration  (ft/sec  ) 

<i>(m)  turntable  angle  at  start  of  mode  m  (deg) 

A  difference  between  bias  determined  in  modes  3  and  4 

and  bias  determined  in  modes  l  and  2  (ft/sec2) 

C.  3  PROGRAM  INPUT 

The  general  format  for  the  input  tape  will  be  as  shown  in  Section 
C.4.  Selected  values  from  each  2-sec  cycle  of  a  60-posit.ion  commutator 
are  recorded  in  digital  form  on  the  input  tape.  Each  commutator  cycle  con¬ 
stitutes  one  data  frame. 

The  signals  of  primary  interest  to  the  data  analysis  program  are: 


a.  Time;  A  nonrepeating  time  signal  from  an  8-day  clock 
aboard  the  satellite  sampled  once  every  2  sec 

b.  Pulse  count:  Signed  accelerometer  pulses  accumulated  for 

1  sec  in  an  up-down  counter  aboard  the  satellite  and  recorded 
twice  every  2  sec 

c.  Mode:  The  experiment  continuously  cycles  through  four  pri¬ 
mary  modes  and  a  transitional  mode.  The  length  of  each  pri¬ 
mary  mode  including  the  appropriate  transitional  mode  is 

512  sec.  The  totnl  calibration  cycle  is  thus  2048  sec 

d.  Pickoff  switch:  The  time  at  which  the  turntable  reaches 
point  A  (forward)  and  B  (aft)  is  recorded  on  the  digital  tape 

e.  Thruster  signal:  A  signal  indicating  which,  if  any,  of  the 
six  attitude  control  thrusters  has  fired  in  the  last  1,2  sec 
is  given  twice  every  2  sec 
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The  calibration  cycle  is  divided  as  shown  in  Table  C.  1.  The  cycle 
is  repeated  every  2048  sec. 


Table  C.  1.  Calibration  Cycle  Modes 


Mode 

Position 

Start  Time, 

sec 

Stop  Time, 
sec 

Elapsed  Time, 
sec 

1 

Low  Speed 

0.  00 

512.  00 

512.  00 

2 

High  Speed 

512.  00 

1024.  00 

512. 00 

5 

High  Speed 

1024. 00 

1031. 38 

7.  38 

3 

Stop  A 

1031. 38 

1536.00 

504.62 

5 

High  Speed 

1536. 00 

1569.  71 

33.  71 

4 

Stop  B 

1569.  71 

2048.  00 

_ 

478.  29 

Time  gaps  in  the  telemetered  data  will  occur  due  to  difficulties  in 
transfering  the  data  recorded  on  board  the  satellite  to  various  ground 
locations. 

Other  inputs  are  the  Mode  1  centripetal  acceleration  (A),  and  the 
ground  calibrated  values  of  scale  factor  (f)  and  bias  (b). 

C.  4  FLOATING  BINARY  FORMAT 

Mode:  Binary 

Density:  800  bits  per  inch 

Record  Length:  Data  records  are  variable  but  are  set  at  117  sets  of 
time  and  quantity  (quantities)  per  record 

Indicative  records  may  contain  up  to  19  three-word 
indicatives 
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rest  aLStdr/tC°rd  °f  f  file  iS  the  indica«v*  record, 
rest  are  data  rec(  is 


All  the 


N^lh<“)irV,Td  °f,a  data  r0C°rd  iS  follow=d  by 
iN  (l  _  N  <  19)  quantity  words 

Number  of  words  per  record  =  42  [T  +  N(Q)] 

£te£2;.  with  ^  indica*wa 

quantUy  ^  ^  ^  “SS  than  42  of  time  and 

There  may  be  any  number  of  files  per  tape 
One  (1)  end -of -file  indicates  the  end  of  the  data  channel 
Jahu  tape  COnSecutivo  '-d -of -files  “dilate  the  end  of  the 


INDICATIVE  FORMAT 
Hit 


1st 

Word 


2nd 

Word 


1  -6 


7-12 


13-18 


19-24 


25-30 


31-36 


1  Character 


3rd 

Word 


Flight 


Channel 


Hits 


Pin 


Orbit 


Stat ion 

i  Rl 


:un 


data  word  format 

Bit 

S  I  Sign 


One 

Word 


A  negative  number  will  have  a  "l"  in  the 
have  a  O  '  in  the  sign  t > t . 


1 


35 


Characteristic 


Mantissa 


sign  bit.  A  positive  number  will 
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N'  Number  of 
Indicatives 
(1  *  N  : .19) 


1  Set  of  T  ime 
and  Quantity 


1  st  Digit 
of  Orbit 


One  Record: 

42  Sets  of  Time 
and  Quantity 

Words /Records 
42  l  T  +  N(Q)] 


T  ime 


C.5  BLOCK  A 


The  1-sec  accumulations  from  the  input  digital  tape  are  processed 
in  the  following  manner.  If  there  is  a  thruster  firing  signal  in  the  first  half 
of  the  commutator  cycle,  then  a  value  of  2048  is  substituted  for  the  original 
value  of  the  pulse  count  for  both  the  first  half  of  the  commutator  cycle  and 
the  second  half  of  the  previous  cycle.  If  there  is  a  thruster  firing  signal  in 
the  second  half  of  the  commutator  cycle,  then  2048  is  substituted  for  both 
halves  of  the  commutator  cycle.  This  is  outlined  in  Table  C-2. 


Table  C-2.  Accelerometer  Pulse  Corrections 


Next  a  4-sec  accumulation  of  accelerometer  pulses  is  made  by 
summing  the  pulse  count  for  four  consecutive  1-sec  intervals,  substituting 
an  average  value  for  those  1-sec  counts  whose  magnitude  exceeds  2000.  In 
case  all  of  the  four  exceed  2000,  an  extrapolation  is  made  from  previous 
4- sec  counts. 

If  n  =  number  of  l-sec  counts  where  !<j’'(i  t  k  -  i)\ >2000,  then  for 
n  <  4,  one  finds 

4 

or'  (i)  -  y"(i  4  k  -  2)  for  |<r"  (i  +  k  -  l,  |  <  2000 
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■IT’/vr-p  '■‘p  T 


For  n  =  4,  one  finds 


<r '  (i)  =  a'  (i  *  4>  +  *'  li  -  8) 

The  above  equation  should  also  be  used  to  bridge  short  «20  sec)  breaks  in 
the  data. 

The  printout  from  Block  A  should  consist  of  every  fifteenth  value 

°  a  along  Wlth  the  corresponding  time  and  mode: 


Time,  sec 


Mode 


Accelerometer  Count 


m(i)  cr'(i ) 

There  will  be  one  sample  per  minute  of  experiment  time.  An  asterisk  m 

should  be  used  to  indicate  those  accelerometer  counts  for  which  n  =  4  and 
the  extrapolation  formula  is  used. 

^  DjEV^ELOPMENT  OF  MODE  1  EQUATIONS 

If  it  is  assumed  that  the  drag  during  a  400-sec  period  may  be  repre¬ 
sented  by  a  second-order  polynomial,  then  the  accelerometer  counts 
accumulated  in  20R  sec  is  given  by 


t+10R 


or  =  f 


J"  £a  +  b  +  (d  +  Dt  +  j  Dt^j  cos (u>t  +  0)  1 
t  - 1  OR  J 


where 


W  .  1 

W  =  TTTTC  sec 


C-7 


Integrating  the  above  equation  results  in 


j  (A  +  b)  +  D  +  (yj  tt2  +  £a/3-1  )"7'lco8(wt  +  0)  -  ( i  -  £\ /3  )“  sin  (wt  +  0) 

w  J 


+  Dt  cos  (ajt  +  0)  -  (l  -  Z^)  ~  8in{wt  +  0)  +  j  Dt^  cos(wt  +  0) 


1  *.2 


This  may  be  further  reduced  by  trigonometric  identities  to 


Sf=i(A  +  b)+  D  cos  $  -  ( 1  .  Lssp.  cos  ut 

u> 

-  [d  sin  0  +  (l  -  Jn/3  C°B-*  +  (77ff2  +  )6  Si2"?  j  sin  wt 

Cl) 

+  C08  0-^1-  ^r»s/3~j^  COS  wt  ■  S*n  0  +  \  1  “  ^J'^'-^jt  8^  wt 


1  ••  2  1  ••  2 
+  jD  cos  0  t  cos  ut  -  J  D  sin  0  t  sin  ut 


In  matrix  form, 


*  y  (A  ♦  b»  ♦  yr  co>  <»*t  *,n  1-1  J  (wt)4,  CO«  wt  J  I  0 


-72  *  I  2  w  s/3  *”  *  |2»  *n  .  72  D  cot  0 

-7  2  -7  2  ♦  }2»  JT  (6/w)  (in  0 

0  72  (D/u*(  coi  0 


i  r  i  2  1 

■jr  jttn  wt  wt  cot  wi  j  (wt)  t in  wt  |  0 


.72  .  12*^1  -  I2«>fj  »  72  D  tin  « 


(D/w)  cot  « 


I  (D/w  )  tin  « 


Allowing  t  to  take  on  the  values 
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t(j)  =  10(2j  -  19)R  *  -  J  -  18 

yields 

wt(j)  =  -g-  (2j  -  19) 


When  the  values  of  t(j)  are  substituted  into  the  above  matrix  equation  and  odd 
and  even  components  are  sepr  "d,  the  even  matrix  equation  becomes 


This  may  be  solved  in  a  least-squares  sense  by  first  premultiplying  by  a 
matrix  with  orthogonal  rows 
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Inverting  the  square  m  trices  yields 


Similarly  the  odd  matrix  equation  becomes 


Premultiplying  to  get  the  least- squares  solution  yields 


1 

-2 

-1 

1 

2 

1 

-1 

-49 

4 

41 

31 

-4 

-22 

-13 

-1590 

-3265 

-1570 

-352 

-611 

-364 

1325 

-2 

4 

2549 


-1 

5 

1329 


o(18)  -  or(l ) 
o(17)  -  <r{2) 
o(16)  -  o(3) 
o( 1 5)  -  o(4) 
0(14)  -  o(5) 
o ( 1  3)  -  o(6) 
o(l2)  -  o(7) 
o(ll)  -  o(8) 
o(10)  -  o(9) 


1 

6 

0 

0 

108 

108  -  6tti/T 

• 

167n2  +  6m/3  -  1Q8 

D  sir  0 

0 

IT 

S’ 

0 

0 

1938  i/T 

-1  9381/1  +  1607ir 

(D/w)  cos  $ 

0 

0 

81431  „ 

“T-rJ 

0 

0 

1 

■ 

(D/id2)  sin  $ 

Inverting  the  square  matrices  yields 


D  am  e 

1  16071.-  3086**1/3 

18  0  0 

(ft/w i  coa  * 

s  m: 

0  6lfl  19  38l/l  -  1607# 

0  5814#  0 

(6/u*)  atn  * 

0  0  1938  <Ti 

0  0  92604426**./! 

-t 

■  2 

.1 

1 

2 

K 

•  49 

4 

41 

31 

-4 

-1590 

.  32  69 

.1970 

-352 

-611 

t  -t 

.21  .11 

•  364  1321 


-2  -t 

4  9 

1949  1329 


»(18)  -  /(II 
»(17)  .  0(2) 

»(lt)  -  0(11 
ff(ll)  -  *(4) 
»(I4)  -  »<9) 
»(1 3)  -  <r(6» 

#02)  -  o(?) 
ulltt.  <r(8| 
#00)  -  »(9) 


C-ll 


iili  ni  I  in  -  -nritfiT n  r- - -  &*^^*** j- -s-»~ .*  »=r 


-46 


C.  7  DEVELOPMENT  OF  MODE  2  EQUATIONS 

If  it  is  assumed  that  the  drag  during  a  400- sec  period  may  be 
represented  by  a  second-order  polynomial,  then  the  number  of  accelerometer 
counts  accumulated  in  20R  sec  is  given  by 

/t+10R 

J^4A  +  b  +  |d  +  Dt  +  2'Dt^jcos(Qt  +  p)j  dt 

t-1  OR 


where  w=  (tt/30R)  sec  .  Integrating  the  above  equation  results  in 

(4A  +  b)  +  |n/"3  D  +  +  y  -  co8(ut  +  0) 

-(V3-  j)  “  &in(ut  +  p)  t  -/T  Dt  cos(wt  +  p) 

-(*/3  -  y)  —  sin(ut  +  0)  +  ‘Jl  j  Dt^  cos(wt  +  0) 

This  may  be  further  transformed  by  trigonometric  identities  to 


we  2u 

T  =T 


sin  wt 


^.<4A*b)  +  [/Tdco.  »-f^-y)^vLj  +  (TV”2^'+T-'f3)^£T^  cos“‘ 

L  <*)  « 

-  r./3D3'n»t(s/3.|)°c°J  t  +  (ig  ,2Vr+5-^)£^J“lsi„  ot 
L  u  • 

+  |»/3D  cos  p  -  -  y  cos  wt  ■  D  sin  p  +  («/3  -  y  )  jt  sin  wt 

I  2  j  2 

+  D  cos  Ot  cos  wt  -  v/T y  D  sin  pt  sin  u»t 
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The  matrix  form  is 


C-  13 


This  may  be  solved  in  a  least-squares  sense  by  first  premultiplying  by  a 
matrix  with  orthogonal  rows: 


Inverting  the  square  matrices  yield 


C  - 14 


Similarly,  the  odd  matrix  equation  becomes 


L 


.1(18) 

-  tr  { 1 ) 

1 

17 

289 

9  (17) 

-  9(2) 

0 

-30 

0 

<r  (16) 

-  <r(3) 

>“ 

1  3 

-169 

iT  ( 1 5 ) 

-  <r  (4) 

1 

1 1 

121 

IT  ( 1 4 ) 

-  c  (5) 

1 

=  18 

0 

-18 

0 

<r(l  3} 

-  t(6) 

.  I 

7 

-49 

it  ( 1 2 ) 

-  <r (7) 

1 

5 

25 

<T  ( i  i ) 

-9(8) 

0 

-6 

0 

-t(IO) 

-  9(9) 

-1 

1 

-1 

54  54  -  6«VT  3tt2  +  6tt^3  -  54 

D  nn  0 

0  6w-J3  -6ffr/7  +  2tt2 

(D/w)  cos  0 

0  0  3tt2 

••  2 

(D/u  )  •  in  0 

Premultiplying  to  get  the  least- square s  solution  gives 


i 

5 

20 


0 

-10 

5 


-I 

5 

-10 


0 

-6 

3 


-I 

3 

6 


1 

1 

■  It 


0 

-2 

1 


-I 

I 

13 


9(18)  -  cr{  I  ) 

<r(l7l  -  <r{ 2) 
<t(1  6)  .  cr(3> 
9(15)  -  cr(4) 
9(14)  -  cr(5) 
9(1  3)  .  9(6) 
<r(12»  -  9(7) 
9(11)  -  9(8) 
«r(10)  -  <r(9) 


w 

7 


1 

T 

0 

0 

54 

54  +  6wi/3 

4 

334"2  -  6tfi/3  -  54 

* 

D  sin  0 

- 

0 

210” 

0 

0 

nT5 

(b/w)  COI  0 

0 

0 

1 1 48ff2 

■ 

0 

0 

l 

« 

(B/w2)  tin  0 

Inverting  the  square  matrices  gives 


D  (in  * 


(B/u.1  col  0 


(D/ u  )  sin 


w 


-  3  */T  -  * 
3  V7 
0 


-328*  +  16it*/I 
-t8if<V3  +  54 

54 


18 

0 

0 


0 

1 890* 
0 


0 

0 

6I992"2 


-1 
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C.  8  BLOCK  D  NUMERICAL  COMPUTATION 


The  calculations  of  Block  B  are  performed  whenever  there  is  a 
complete  mode  or  a  beginning  fragment  of  Mode  3  or  4.  The  start  of  each 
mode  is  determined  by  a  change  in  the  mode  signal.  For  Modes  3  and  4, 
the  starting  time  is  considered  to  be  the  start  of  the  previous  transitional 
mode  (Mode  5).  The  starting  time  can  be  determined  either  by  the  change 
from  the  previous  mode  to  the  transitional  mode,  or  by  the  change  from 
Mode  5  to  Mode  3  or  4  minus  a  known  correction.  The  time  corrections  are 
outlined  in  Table  C-3. 


Table  C-3.  Mode  Time  Correction 


The  starting  time  (t^)  is  determined  by  subtracting  the  time  correction  from 
the  time  a  change  in  mode  is  detected.  If  a  determination  has  been  made  at 
the  start  of  Mode  5,  the  redundant  determination  at  the  end  of  Mode  5  should 
be  omitted. 

For  Modes  /  and  2,  the  4-sec  accumulations  of  accelerometer 
counts  are  combined  into  18  accumulations  in  the  20R  -  sec  intervals  from 
+  300  -  180R  to  tg  +  300  r  180R;  therefore. 


C  - 16 


17 

<r(j)  =  ^  g(k)  <r'(t0  +  74  +  2 0 j  +  4k)  1  £  j  £  18 

k=  1 

whe  re 

g(k)  =  1  for  5  ( j  -  9  )R  +  56  2  5j  +  k  2  5(j-10)R  +  57 

g(k)  =  5{j  -9)R  +  57  -  5j  -  k  for  5(j-9)R  f  57  >  5j  +  k  >  5 ( j -9)R  +  56 

g{k)  =  5j  +  k  -  56  -  5{j-10)R  for  5  ( j  -  1 0)R  +  57  >  5j  +  k  >  5(j  - 1 0)R  +  56 

g(k)  =  0  for  5(j-10)R  +  56  2  5j  +  k 

5 (j  -  9)R  +  57  s  5j  +  k 

Numerical  calculations  can  best  be  performed  by  first  forming 
auxiliary  summations  as  indicated  by  the  following  matrix  equations.  For 
data  from  Mode  1 


C  - 1 7 


C  - 18 


'  SS(2,  1)  r  18  0  0 

SS(2,  3)  =  0  1  890tt  0 

.  SS(2, 5)  J  0  0  61  992it2 


‘  1  0-110-1  1 
x  5  -10  5  3  -6  3  1 

20  5  -10  0  3  6  -11 


If  data  are  available  to  compute  the  summations  for  consecutive  Modes  1 
and  2,  then  calibration  calculation;-  are  performed  to  update  the  values  for 
accelerometer  scale  factor  (f)  and  bias  (b) 

f  =  TMra{2SS(2,0)  ’  SSf1’0*  -  2n/I[SS(2,4)  +  3SS(1,4)] 

+  ttn/T[2SS(2,  6)  -  35SS(1, 6)]| 

b  =  90lj{2SS(1’0)  -  SS(2’  °>  +  ^  [l2SS(l ,  4)  +  SS(2,  4)] 

+  vJJ  [70SS(i,  6)  -  SS(2,  6)]} 

If  either  Mode  1  or  2  is  not  complete,  then  the  previous  values  of  scale 
factor  and  bias  are  used  to  compute  drag  acceleration  and  phase  angle. 


For  Mode  1 


0(1)  =  f[SS(l.l)  -  (6  VT-  rr)SS(l,  3)  -  tt(3086tn/T-  1607)SS(1,  5)j 


ji  +  [JSiU 
(  L  ss  ( i ,  i 


2)  -  (5tt  -  6nT3)SS(1,4)  -  tt(886tt  -  35nT3)SS(1 
)  -  (6 n/T-  tt)SS(  1 ,  3)  -  *(3086 tfs/T-  l607)SS( 


(1 .  5)  J  J 


6(1)  =(^)  f  l  6SS (1.4)  -  (19W3  -  3  ir)SS(  1 , 6)J 

v  |  3  +  ~6n/3SS(1,  3)  -  (1  607tt  -  1  938>/T)SS(l,  5)' 2  j1 
(  L  6SS(1, 4)  -  (194*/!-  35n)SS(l,  6)  J) 


p,,  v  _  /  "  \J  1938  ,  L  97SS(1 , 6)  ' 

D()"l30R/  ~T~  SS(  ’  5)  j  3  (.32  3SS  ( 1 ,  5) . 


p(  1 )  -  arc  t 


For  Mode  2 


an  ~SS(1-  U  ~  ff)SS(1»  3)  -  ff(3086TTNT3  -  1607)SS(1,  5)  ~j 

L  SS (1,2)  -  (5*  -  6\H>  )SS(1 , 4)  -  *(886*  -  35n/T)SS(1  ,  6)  J 


0(2)  =  -•[  SS(2,  1)  -  (3n/T+  *)SS(2,  3)  -  2*(i64*  .  9'/3)SS(2,  5)J 


X  \  1  + 


SS(2,  2)  -  (tt  .  yT)SS(2,  4)  -  *(52*  4  nT?)SS(2,  6)  1 2 11  ! 

„SS(2,  1)  -  (3n/T  +  *)SS(2,  3)  -  2tt(  1  64tt  .  9\f3)SS(2,  5)  J  J 


6(2)  -  (^~)  ^lSS(2.4)  -  (3vT  4  *)SS(2,  6)] 


(h9^ 

f  LsS(2, 


3)  -  6 (r  -  \TT )SS(2 

4)  -  (  3n/T  4  *)SS(2 


llif!1 

,6)J  ) 
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6“>=(W)  ?SS(2.S>)36+[iJ|t|»J  f 

*(2)  =  arc  tan  pS^.1)  -  (3>/T  +  *)SS(2,  3)  -  2n(\64n  -  9\TT)SS(2.5) 
L  SS(2,  2)  -  (tt  .  n/T )SS(2,  4)  -  Tr(52Tr  +  -/T)SS{2,  6) 


For  both  Modes  1  and  2 


d(j)  =  D  +  1 0 (2 j  -  21)D  +  50(2j  -  21)^6  1  <  j  <  20 


For  both  Modes  3  and  4 


(J)  =  Z  tT'  (t0  +  78  +  20j  +  4k)  1  <  j  <  20 


For  Mode  3 


d(j)  =  loir  *  b 


»<3>  =  *b  +  jrsrhom  Z  {738, 507  ' 24> 519j  +  i3’ 663j2),r{j) 

j=i 


For  Mode  4 


d(j)  =  b  - 


0(4)  =  b  - 


(-1  68,  245  +  85,  349j  -  4,  257j2)cr(j) 
J--1 
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For  Modes  1  ?  -i  ,  . 

'  '  '  an  .the  time  of  each  drag  value  is 

=  tQ  f  90  +  20 j 

The  values  of  d(i)  should  k 

Whenever  consecutive  Modes  3  and  \  C0™PUted  f°r  Part«l  Mode  3  or  4. 
he  computed  as  follows  complete,  the  bias  check  A  should 


A  =  \  fD(3)  -  D(4)] 

A  printout  should  be  made  for  earh  f 

-ale  factor.  bias.  phase  and  ^  ^  * 

Of  these  quantities  will  be  available  Th  f  ’  "  “C°n’plete  cycles  not 

as  shown  below;  '  6  forrnat  of  the  printout  shall  be 


Time, 


tg  (Mode  1 )  or 
tQ  (Mode  3)  -  1024 


Scale  F actor 

P^£/ft/se^  ft/sec2  VY/s^c^’  *'  Phase  2- 

- - -  - ~ dec 


In  addition.  a  printout  of  the  20  vai,  ,  , 
made  for  each  mode.  Less  than  20  valu  ^  °  rag  acceleration  should  be 
MOd*  3  «  "  Th.  format  is  ^ 


I  irne,  sec 


ft /sec' 


,J  d  (j )  /  32.  145 

The  values  of  time  and  acceleration  u 

can  be  used  as  the  ^  *  --q:::0:;:er;::i:;d  in  a  form 
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Block 

B 


^  Block  B 
Printout 


C  -23 


APPENDIX  D 

LOGACS  DATA  ANALYSIS  PROGRAM 


NOTE:  A  complete  development  of  the  computer  program  and  subroutines 
used  for  the  LOGACS  data  analysis  is  presented. 
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appendix  e 


1 


LOGACS  SENSED  ACCELERATION  ALGORITHM 


The  LOGACS  accelerometer's  proof  mass  is  connected  to  the  Agena 
ve  tele  through  a  force  suspension  system  in  which  one  component  (defined 
as  the  input  axis  component,  is  measured  explicitly.  This  appendix  develops 
the  algorithm  for  this  measured  acceleration  component. 


Assuming  the  Agena  vehicle  is  a  rigid  body,  a  force  vector  diagram 

of  the  vehicle/accelerometer  system  in  a  fixed  reference  frame  can  be  con- 
structed  as  shown  in  Fig.  E-i, 


f  ig.  E-I.  Force  Vector  Diagram 


The  symbols  in  Fig.  E-l  are  defined  as  follows: 


R.  =  radius  vector  from  the  earth  center  to  i  (i  =  v,  a) 

=  vector  sum  of  external  forces  applied  to  i  (i  =  v,  a; 

F  represents  accelerometer  bias  term) 

^  d 

Fgi  =  suspension  force  vector  (i  =  v,  a) 
m.g(R.)  =  gravity  force  vector  (i  =  v,  a) 

From  these  definitions,  the  acceleration  equation  for  each  mass  is 

K  ~~  +  I(R v)  (E-l) 


+  F 

sa 


m 


a 


+  *<Ra> 


Using  the  relationships 


(E-2) 


1  =  R  -  R 
a  v 


R  =  R  -  R 
a  v 


F  +  F  =0 
sv  sa 


and  substracting  Eq.  (E-l)  from  (E-2)  results  in 


..  F  +  F  F  +  F 

rT  _  a  sa  v  sa  ,  — /sr  .  — , 

R  = - — - — - +  g(R  -  g  R 

m  m  °  a  °  v 

a  v 


(E-3) 


.  . -  -  -  --  --i..  ; — - — - -II  lltfililtiJllitlill'  fnilif 


Assuming  the  accelerometer  mass  is  negligible  with  respect  to  the 
vehicle  mass,  the  acceleration  sense*!  by  the  accelerometer's  force  suspen¬ 
sion  system  is 


=sr  +  «-[8<R  a»  -  «(RV>]  "sr  (E-4) 

a  v  a 

The  acceleration  vector  R  expressed  in  the  rotating  vehicle  axis  frame  is 

R  =  r  +  ux(uxr)+uXr  +  2uxr  (E- 5) 

« 

where  w  and  ware  the  inertial  angular  rate  and  acceleration  vectors, 
respectively,  of  the  vehicle  axes  and  r  is  the  displacement  vector  expressed 
in  body  axes  coordinates.  Neglecting  the  bias  term  (F  /m  ),  changing  the 
symbols  of  the  force  to  mass  terms,  and  using  the  symbol  g(r)  for  the  gravity 
gradient  term  results  in  the  following  equation  for  the  sensed  acceleration 
vector  expressed  in  body  axis  coordinates 

•  • 

A  =  A  +  7  + u  X  (7  X  7)  +  w  X  7  +  2u  X  r  -  g(r)  (E-6) 

a  eg 


where 


r  = 


X  +  r  cos  0, 
ag  a  rt 


Y  +  r  sin  6 
ag  a  rt 


i  z 
\  a* 


instantaneous  coordinates  of  proof 
mass  with  respect  to  the  vehicle  eg 
(see  Fig.  19) 


u 


u> 


w 


w 


Inertial  body  rates 


g(r)  =  gravity  gradient  term  (see  below) 


E-3 


. . 


Since  the  coordinates  of  the  LOGACS  package  with  respect  to  the 


eg  (X  ,  Y  ,  Z  )  and  the  turntable  radius  arm  r  are  assumed  constant, 
ag  ag  ag  a 


the  relative  velocity  and  acceleration  of  the  proof  mass  are 


-sin 

r  -  r  cos  <p. 
at  1 


(E-7) 


cos  0t\ 


'-sin  0t' 


-  ra*t  sin  0t  +  ra0t  cos  0t 


(E-8) 


(E-9) 


The  gravity  gradient  term  is  developed  assuming  a  spherical  earth 
as  follows: 


i(R  )  =  -  R 


3  v 


g(6R} 

RJ 

v 


—  —  t 

3R  R1  , 

— - 1  SK  -  j  0 

R“  ° 

v  ■*  0 


-i  0  UR 


3 


where 

0,v|/,V>  =  ordered  rotations  (pitch,  yaw,  roll)  of  the  Agena  vehicle  axis 
with- respect  to  the  orbit  plane  axis. 

When  it  is  assumed  that  the  angles  (0,i}/,v>)  are  small  (cos  =  1, 
sin  =  angle,  sin  products  =  0),  the  order  of  rotations  is  immaterial  and  the 
transformation  reduces  to  a  skew  symmetric  matrix 


In  the  orbit  plane  coordinate  system,  the  z  axis  is  directed  along  the  negative 
radius  vector,  the  y  axis  is  directed  normal  to  the  orbit  plane,  and  the  x  axis 
is  directed  forward,  formating  a  right-hand  system. 


E-5 


.  r-.r-r  -  I,  frifl  .  mri 


f-w*^pwr»Tr-«  --^>-  <•->  -  ■''—••  s-;  ’  -rrrv-^  S^V^^  ' ‘^  ^mv  ”*1 


With  the  above  definitions,  the  gravity  gradient  term  ij 


g(f)  =  Mbog(6R) 


-i  o  o 


=  wKo  0  -1  °IM 


■  .7**.. 

0  0  2 


T  — 

Mio* 


-i  o  -3e" 


=  Uq I  0  -i  3^1 r 


-3e  3<fi  2 


(E-  10) 


Substituting  Eqs.  (E-7),  (E-8),  (E~9),  and  (E- 10)  into  Eq.  (E-6) 
and  collecting  terms  results  in  the  following  x  and  y  components 


Asx  =  Axcg  “  ra^cos  *t^t  +  Wz*  +  wy  "  W0J  +8in  +  "z  "  wxwy)| 


2,2  2 


“Xag<wy  +wz  *  uo>  “  Yag<“*  ‘  wx<V  +  Zag(Vz  +  «y  +  ^o9* 


Lsy  =  Aycg  "  rJsin  *t^t  +  wz)2  +  wx  "  wo]  +  cos  +t{~*t  “  “z  ’  WXM 


-Xag<“x“y  -  “z>  -  Yag'“x  +  “a  *  “«•  *  Zag<“x  *  “y“*  +  3“o*’> 


The  input  axis  acceleration  to  the  LOGACS  accelerometer 


Asi  =  Asx  COS  *t  +  Asy  *in  *t 


Performing  the  above  transformation  and  collecting  terms  results  in 

Asi  Axcg  cos  +  Aycg  sin  "  raC(*t  +  w z ^  +  (cos  +tUy  "  sin^twx)2  * 
+  Xag  C8in  +  “xV  -  cos  *t<wy  +  “l  " 

+  Yag  t-COS  0t{wxwy  "  “z>  "  sin  *t(wx  +  uz  “  wo^ 

-f*w.  r 

+  Zag[cos  *t(«x<-z  +  +  3^6)  -  sin  ^  -  ty^  +  3u>^)] 


(E-ll) 
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Corrected  Raw  Accelerometer  Data 
Revolution  54  Perigee) 


Fiy.  F-4.  Corrected  Raw  Accelerometer  Data  (Modes  1  and  Ly 
Revolution  54  Perigee) 
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APPENDIX  G 


REDUCED  D>'r A  PLOTS 


This  appendix  present*-  .he  plots  of  reduced  data  from  the  LOGACS 
flight.  Figure  G~1  presents  the  drag  acceleration  profile  reduced  from  the 
accelerometer  data  and  the  drag  acceleration  profile  from  the  best-fit 
ephemeris.  Also  included  are  plots  of  latitude,  longitude,  and  altitude  from 
the  besi-fit  ephemeris.  The  orb.t  rev.  number  is  indicated  .'.ear  each 
perigee.  Figure  G-l  is  separated  into  four  plots  one  for  each  lay  of  the 
mission  [Fig.  G-l(a)  is  for  23  May,  Fig.  G-l(o!  for  24  May,  et  .  ]. 

Figure  G-2  presents  information  relating  to  the  contro’  system. 

The  drag  acceleration  profiles  rc  ilso  ircluded  sc  that  correla  ions  of 
anomalous  responses  can  be  observed.  Figure  G-2  is  also  separated  into 
four  plots  and  the  rev.  number  is  indicated  near  perigee. 
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(a)  Acceleration  Profile  and  Ephemeris  Data  (May  £3) 
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